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NOTATION 


A  dummy  integration  variable  for  area 

Aci  surface  area  between  the  mass  of  material  compacted  by  the  landing 

surface  and  the  mass  of  uncrushed  crushable  material 

Acih  horizontal  planar  projection  of  Aqi 

^po  surface  area  of  lower  portion  of  payload  over  which  shear  and/or 

normal  stresses  act  to  decelerate  payload 

^poh  horizontal  planar  projection  of  Ap^  and,  in  the  absence  of  shear 

effects,  Api 

Api  surface  area  between  the  mass  of  material  compacted  by  the  payload 

and  the  mass  of  uncrushed  crushable  material 

Agj^  horizontal  cross-sectional  area  in  the  volume  of  figures  3 

and  4  (the  volume  that  would  lie  beneath  the  landing  surface  if 
there  had  been  no  crushing) 

b  dimensional  dummy  integration  variable  used  in  equation  (A27) : 

E  /r2  -  r2 

bp  dimensional  dummy  integration  variable  used  in  equation  (A29) ; 

E  /Rp2  -  rp2 

e  dimensionless  variable  describing  the  relative  motion  of  payload 

•ip  ■  R 

penetration;  defined  as  — -  in  equation  (A39) 

P 

©max  maximum  value  of  e  during  payload  penetration 

Rp 

Pcl(y)  identical  to  Fd  (z)  if  -g— y  is  substituted  for  z 

Pci(^)  dimensionless  vertical  crushing  force  acting  over  area  and 

tending  to  decelerate  the  mass  above  ;  defined  in  equation 
(A36) 

dimensionless  vertical  crushing  force  acting  over  area  Ap^  and 
tending  to  decelerate  the  mass  above  Ap^;  defined  in  equation 
(A45) 

g  general  symbol  for  acceleration  due  to  gravity 

ge  value  of  g  on  earth  (32.17  ft/sec^  herein) 

gp  value  of  g  at  the  landing  site  (12.3  ft/sec^  herein) 


V 


local  height  of  the  volume  which  would  lie  beneath  the  landing  sur¬ 
face  if  there  had  been  no  crushing;  the  local  height  of  V^;  see 
figures  2,  3,  and  4 

local  height  of  the  compacted  volume  of  material  compacted  by  the 
landing  surface;  the  local  height  of  mci;  defined  in  equation 
(A12);  see  figures  2,  3,  and  4 


local  height  of  the  compacted  volume  of  material  compacted  by  the 
payload;  the  local  height  of  mpi;  defined  in  equation  (A12) ;  see 
figures  2,  3,  and  4 

/getSEA)\/EdW  \ 

defined  in  equation  (17)  as  - ^ —  —  I - 

\  Uo"  /V'"/VV  "“/ 

dummy  integration  variable  for  z  in  equations  (A3 7)  and  (A49) 


defined  in  equation  (A39)  as 


TrRp2ao 

^o%dSe 


defined  in  equation  (A33)  as 


^o^md^( 


distance  between  the  fictitious  compacted  region  determined  by  ej 
on  the  landing  surface  and  that  on  the  payload  (or  the  payload 
itself,  in  the  absence  of  penetration);  evaluated  in  equation 
(5) 

variable  used  as  di  to  define  a  length  element  of  a  rod  of 
crushable  material;  used  in  developing  equation  (All) 

original  mass  of  uncrushed  crushable  material;  see  figures  2(a), 
3(a)  ,  and  4  (a) 

mass  of  material  compacted  by  the  landing  surface;  mass  of  ; 
see  figures  2  through  4 

mass  of  payload;  see  figures  2  through  4 

mass  of  material  compacted  by  the  payload;  see  figures  2(c),  3(c), 
4(c) 


first  time  derivative  of  mi 


dimensionless  mass  defined  in  equation  (A34)  as 


^  "^co  '  ^cl 


m  (y,e)  dimensionless  mass  defined  in  equation  (A44)  as 


^co  “  ^cl 


VI 


dimensionless  mass  defined  in  equation  (A42)  as 


m  (e) 
pen 


/  ^0 

defined  in  equation  (23)  as  —  I  1  +  :r - - ^ 


defined  in  equation  (24)  as 


1  +  (m^o/mpo) 

[l  +  CUo^/2ednp  geRp)] 


defined  in  equation  (25)  as  -  (see  also  eq.  (26) 

for  payload  penetration  with  the  simplified  design) 


desired  value  of  n 


Pmax 


^max^^Pmax 


maximum  permissible  value  of  (2000  herein) 

g  loading  with  Upg^  as  the  payload  deceleration 

maximum  value  of  np  that  occurs  during  impact  stroke 

absolute  displacement  of  uncrushed  crushable  material  after 
initial  contact  with  landing  surface;  see  figures  2  through  4 

first  time  derivatives  of  q  and  qp 

second  time  derivatives  of  q  and  qp 

absolute  displacement  of  payload  after  initial  contact  between 
crushable  material  and  landing  surface;  qp  =  q  in  absence 
of  payload  penetration;  see  figures  2  through  4 

values  of  q  and  qp  at  end  of  impact  stroke 


value  of  q  (q  =  qp)  at  which  payload  penetration  occurs  or 
would  occur  if  unbonded 

overall  radius  for  spherical  system 

payload  radius  for  spherical  system 

polar  radial  coordinate  of  point  in  A^i,  where  coordinate  is 
measured  in  a  horizontal  cross  section  for  spherical 
geometry;  see  figure  4(b)  and  equation  (A25) 

polar  radial  coordinate  of  point  in  Ap;^ ,  where  coordinate  is 
measured  in  a  horizontal  cross  section  for  spherical 
geometry;  see  figure  4(c) 


SEA 


V 


^max 


cl 


.Vci 


max 


V. 


V 


pi 


V, 


specific  energy  absorption;  energy  absorbed  per  unit  weight  of 
the  absorber;  defined  in  equation  (8)  in  terms  of  the 
variables  used  herein 

dimensionless  dummy  integration  variable  used  in  equation  CA36) ; 
^  b 
”  R 

dimensionless  dummy  integration  variable  used  in  equation  (A45) ; 
-  Rp 

time;  used  for  time  derivative  interchangeably  with  dot 
velocity  of  the  uncrushed  crushable  material;  defined  in 
equation  (A16)  as  q 


dt 


value  of  U  at  the  instant  the  crushable  material  hits  the 
landing  surface  (300  ft/sec  herein) 

volume  that  would  lie  beneath  the  landing  surface  if  there  had 
been  no  crushing;  shown  in  figures  3  and  4;  defined  in 
equation  (A14) 

maximum  values  of  and  Vqi,  that  is,  the  values  reached  at 
the  end  of  the  impact  stroke 

compacted  volume  of  material  compacted  by  the  landing  surface; 
the  volume  of  m^^^ 

volume  swept  out  by  the  payload  as  a  result  of  the  relative 
motion  of  penetration;  defined  in  equation  (A14) 

compacted  volume  of  material  compacted  by  the  payload;  the 
volume  of  mpi 

dimensionless  velocity  of  the  uncrushed  crushable  material; 


defined  in  equation  (A54)  as 


dq/dt 


U 


^%dSe^p  ^^md^e^p 


dx 


w 


value  of  V  at  the  instant  the  crushable  material  hits  the 
landing  surface 

dimensionless  velocity  of  the  payload;  defined  in  equa- 


^Ary,^  _  de  dy  _ 

tion  (A54)  as  —  .  .  =  ev  -r:—  +  = 

^^mdSe^p 


dy. 


dx  '  dx  dx 


value  of  V  at  the  start  of  payload  penetration,  if  any 


viii 


original  earth  weight  of  uncrushed  crushable  material,  m^oge 
earth  weight  of  payload,  nip^g^ 

dimensionless  velocity  of  the  uncrushed  crushable  material; 


defined  in  equation  (A33)  as 


^mdSe^ 


value  of  w  at  the  instant  the  crushable  material  hits  the 
landing  surface 

value  of  w  at  the  start  of  payload  penetration,  if  any 
dimensionless  time;  defined  in  equation  (A39)  as  t  / 


^max^^Pmax 


dimensionless  displacement  defined  in  equation  (A39)  as 
Rp  "  ^  Rp 

dimensionless  payload  displacement  defined  in  equation  (A54)  as 

qp  _ 

^  E  e^e  +  y 

values  of  y  and  yp  at  end  of  impact  stroke 
value  of  y  at  the  start  of  payload  penetration,  if  any 
dimensionless  displacement  defined  in  equation  (A33) 
value  of  z  at  np  = 

value  of  z  at  end  of  impact  stroke 


value  of  z  at  the  start  of  payload  penetration,  if  any 

angle  between  normal  to  stressed  area  and  direction  of  maximum 
normal  crushing  stress 

^Pmax 


compacting  strain  of  crushable  material;  more  detailed 
definition  given  following  equation  (1) 


IX 


Pck 


P  cm 

Pp 


PpR 


a 


^po 

^pok 

^v' 

^vpo 

p 


fictitious  value  of  assumed  for  design  purposes;  always  less 

than  0^ 

compacting  strain  when  it  is  uniform  throughout  hypothetical 

crushable  material;  determines  the  surface  for  stress  evaluation 
and  can  be  specialized  in  governing  equations  for  simplifying 
approximations 

value  of  eic  for  actual  crushable  material  (0.8  herein) 

angle  shown  in  figure  4(c)  and  used  in  derivation  of  equation  (A28) 

dummy  integration  variable  for  y 

uniform  density  of  hypothetical  crushable  material  in  swept-out 
volumes  such  as  Vc  in  figure  4;  Pck  can  be  specialized  as  zero 
to  neglect  variable  mass  without  implying  a  massless  crushable 
casing 

uniform  density  of  actual  crushable  material 

payload  packaging  density;  defined  in  equation  (A43)  as 
”^0 

(4/3) 

defined  in  equation  (A35)  as  - 

(4/3)TrR3 

"mostly  normal"  and  "mostly  static"  crushing  stress  of  crushable 

material,  averaged  over  maximum  possible  stroke  prior  to  compacting 

maximum  of  a  as  defined  in  equations  (2)  and  (A24) 

normal  stress  on  the  payload  prior  to  penetration 

OpQ  under  assumption  Opo  =  ^-pok  ~  constant 

vertical  component  of  o 

dynamic  value  of  as  defined  following  equation  (A7) 

vertical  component  of  normal  stress  on  the  payload 
angle  shown  in  figure  4(b)  and  used  in  derivation  of  equation  (A27) 
angle  shown  in  figure  4(c)  and  used  in  derivation  of  equation  (A28) 
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APPROXIMATE  ANALYTICAL  MODELS  FOR  LANDING  ENERGY  ABSORPTION, 
INCLUDING  THE  EFFECT  OF  PENETRATION  BY  THE  PAYLOAD 
INTO  ITS  CRUSHABLE  CASING 
Robert  W,  Warner 
Ames  Research  Center 


SUMMARY 


Two  approximate  analytical  models  are  defined  for  a  landing  configuration 
in  which  a  spherical  payload  can  sometimes  penetrate  into  its  crushable  casing. 
Results  for  both  models  are  found  to  agree  reasonably  well  with  two  previous 
experimental  measurements.  Design  examples  are  presented  for  an  impact  veloc¬ 
ity  of  300  ft/sec.  These  are  based  on  choices  of  zero  or  "perfect"  payload 
bonding,  and  of  either  a  balsa-like  or  honeycomb -like  class  of  crushable  mate¬ 
rial.  The  greatest  difference  between  the  models  for  these  examples  is  a  29- 
percent  discrepancy  in  the  required  maximum  crushing  stress.  A  particular 
pair  of  examples  gives  the  unexpected  result  that  penetration  can  provide  a 
decrease  in  crushable  material  weight  by  a  factor  greater  than  4  when  the 
honeycomb -like  class  of  material  is  required  without  penetration,  but  the  more 
efficient  balsa-like  class  is  feasible  with  penetration. 


INTRODUCTION 


One  means  proposed  for  providing  information  on  lunar  and  planetary 
surfaces  consists  of  an  unmanned  instrumentation  system  that  is  hard  landed 
(with  or  without  terminal  guidance)  but  designed  to  survive  the  impact.  For 
such  a  landing,  a  crushable  casing  is  one  means  for  absorbing  landing  energy 
so  the  payload  can  survive  and  transmit  information  during  and  after  the 
impact.  The  advantage  of  this  approach  depends  greatly  on  how  light  the 
crushable  casing  can  be  made  for  a  given  impact  velocity. 

Various  aspects  of  the  design  of  crushable  casings  have  been 
investigated  in  references  1  through  9,  but  without  including  the  effect  of 
penetration  by  the  payload  into  the  casing.  The  primary  purpose  of  the  pres¬ 
ent  paper  is  to  evaluate  that  effect  analytically.  For  this  purpose,  fairly 
general  equations  of  motion  are  derived  for  the  payload  and  the  crushable 
material.  These  equations  are  then  specialized  for  two  approximate  analytical 
models  in  which  the  payload  and  casing  are  spherical.  The  analytical  models 
will  be  used  in  a  number  of  design  examples  so  that  the  effect  of  payload 
penetration  can  be  evaluated  and  also  one  model  compared  with  the  other.  The 
two  models  will  also  be  compared  with  earlier  models  that  do  not  permit 
penetration  and  with  the  test  results  of  reference  10  and  a  private 

communication.  ^ _ 

^Donald  R.  Cundall,  December  1967.  ' 


OUTLINE  OF  THEORY 


Properties  of  Typical  Crushable  Material 

Material  properties  are  an  essential  input  to  the  theoretical 
development.  A  typical  crushable  material  for  landing  impact  energy  absorp¬ 
tion  has  a  stress-strain  curve  similar  to  curve  ABODE  in  figure  1(a).  The 
material  is  elastic  from  point  A  to  point  B,  and  there  is  rebound  between 
points  D  and  E.  (The  rebotmd  should  be  as  small  as  possible.)  Between  points 
B  and  D  there  is  a  large  volume  change;  and  the  crushing  stress  is  relatively 
constant,  with  an  average  value  a  indicated  by  the  horizontal  dotted  line 
in  figure  1(a).  Point  C,  where  the  crushing  stress  begins  to  rise  abruptly, 
is  called  the  compacting  strain  e  and  ranges  from  0.6  for  close-packed 
crushable  materials  to  nearly  1.0  for  open  crushable  structures. 

The  area  enclosed  by  the  stress-strain  curve  ABCDE  is  the  energy 
absorbed  per  imit  volume  of  crushable  material.  This  energy  is  maximized  for 
a  maximum  permissible  crushing  stress  (which  determines  the  maximum  landing 
vehicle  deceleration  for  a  given  configuration)  if  the  stress-strain  curve 
approaches  a  rectangle.  Hence  the  material  is  not  strained  beyond  point  D  in 
figure  1(a),  where  the  stress  is  equal  to  the  prior  maximum  (point  B) ,  even 
though  the  stress  could  go  higher  as  shown.  If  the  initial  peak  at  B  is  too 
high  to  approximate  a  rectangle,  it  can  sometimes  be  reduced  by  precrushing. 

A  less  desirable  alternative  is  to  accept  the  stress-strain  curve  but  modify 
the  load-deflection  curve  by  changing  the  shape  of  the  crushable  material. 

Stress-strain  and  load-deflection  curves  for  balsa  wood,  plastic  foam, 
and  honeycomb  are  given  in  references  11  through  14.  These  curves  have  been 
determined  by  dynamic  crushing  tests  and  by  nearly  static  tests  of  specimens 
having  uniform  cross  sections.  Figure  1(b)  is  a  sketch  of  a  crushing  test  in 
which  the  specimen  is  compressed  uniformly  across  the  cross  section  by  a 
plate,  and  figure  1(c)  shows  a  test  in  which  the  specimen  is  penetrated  by  a 
plunger.  In  both  tests,  the  material  crushes  in  layers  at  a  loaded  surface, 
which  may  be  at  either  end  of  the  specimen  in  the  case  of  the  plate  loading. 

It  should  be  noted  that  a  certain  (exaggerated)  amount  of  material  is 
shown  trailing  outboard  of  the  plunger  in  figure  1(c)  due  to  shear  resistance. 
This  shear  effect,  as  well  as  friction,  causes  a  difference  in  the  crushing 
loads  of  figures  1(b)  and  1(c).  In  the  case  of  balsa  wood,  however,  refer¬ 
ence  7  indicates  that  this  difference  is  small,  while  reference  13  indicates 
that  it  ranges  from  5  to  15  percent  (for  a  range  of  plunger  sizes). 

It  should  also  be  noted  that  dynamic  and  nearly  static  crushing  tests 
give  different  results.  In  fact,  references  11  and  13  suggest  ratios  of 
static  to  dynamic  crushing  stress  from  0,69  to  0.73  for  various  materials. 
Since  the  maximum  velocity  in  the  dynamic  tests  was  108  ft/sec,  the  dynamic 
effect  is  probably  not  due  to  variable  mass  (i.e.,  the  accumulation  of  crushed 
material  at  a  loaded  surface)  but  rather  due  to  damping  and  dynamic  buckling 
(i.e,,  coupling  between  vertical  and  horizontal  velocity).  Unfortunately, 
the  effects  of  higher  velocities  on  damping  and  dynamic  buckling  are  not 
established  for  the  present  materials. 
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Summary  of  Fundamental  Assumptions  and  Limitations 

In  contrast  to  the  typical  material  just  discussed,  the  material 
assumed  for  theoretical  analysis  has  a  perfectly  rectangular  stress-strain 
curve  (a  so-called  ”rigid  plastic"  shape) .  This  curve  is  represented  by  the 
dashed  lines  in  figure  1(a).  It  is  shown  bounded  by  the  compacting  strain  e 
and  the  average  crushing  stress  o  of  the  typical  material.  In  this  case, 
the  energy  absorbed  per  unit  volume  will  not  match  that  of  the  typical  material 
perfectly,  but  the  boundaries  of  the  rectangle  can  be  adjusted  slightly  if 
need  be. 

There  are  a  number  of  other  analytical  assumptions  that  pertain  to  the 
crushable  material,  and  there  are  several  that  do  not.  For  convenience,  all 
fundamental  assumptions  and  limitations  are  listed  in  the  present  section  as 
follows : 

1.  Rebound  is  assumed  absent, 

2.  It  is  assumed  that  the  effects  of  shear  resistance,  end  fixity,  and 
Poisson *s  ratio  are  adequately  incorporated  in  the  analysis  because 
of  their  presence  in  the  crushing  tests  that  determine  the  so-called 
"mostly  normal"  and  "mostly  static"  crushing  stress  (a,  a^,  and  o^; 
see  Notation) . 

3.  The  "mostly  normal"  and  "mostly  static"  crushing  stress  a  is  also 
assumed  to  incorporate  the  dynamic  effects  of  damping  and  dynamic 
buckling.  This  means  that  dynamic  tests  should  be  used  to  determine 
o  (or  static  dynamic  ratios,  such  as  those  deduced  earlier  from  ref¬ 
erences  11  and  13,  should  be  applied) .  The  velocities  in  the 
dynamic  tests  should  not  be  large  enough,  however,  to  cause  signifi¬ 
cant  variable-mass  effects,  which  are  incorporated  (when  desired)  in 
the  present  equations  of  motion. 

4.  Shear  deformations,  such  as  the  trailing  of  material  outboard  of  the 
plunger  in  figure  1(c),  are  neglected. 

5.  Separate  vertical  rods  of  material  are  assumed  to  crush  vertically 
in  the  energy-absorbing  process. 

6.  The  material  is  assumed  to  compress  to  the  same  compacted  strain 
along  any  axis,  regardless  of  the  axis  of  maximum  normal  crushing 
stress . 

7.  It  is  assumed  that  each  particle  in  the  uncrushed  crushable  material 
is  moving  at  the  same  vertical  velocity  at  a  given  instant  and  that 
the  crushed  material  also  has  a  uniform  particle  velocity,  but  a  dif¬ 
ferent  one  from  that  of  the  uncrushed  material.  This  implies  the 
following  corollary  assumptions: 

a.  Each  successive  layer  of  crushable  material  undergoes  a 
jump  in  velocity  as  it  moves  from  the  uncrushed  to  the 
crushed  region. 
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b.  The  elastic  stress  waves  that  established  the  uniform 
velocities  must  travel  far  faster  than  those  velocities 
(i.e.,  the  uniform  velocities  are  subsonic). 

c.  The  deformation  waves  resulting  from  the  elastic  stress 
waves  must  be  small  enough  not  to  affect  the  uniform 
velocities . 

8.  For  landing  geometries  in  which  there  is  doubt  as  to  where  the  crush¬ 
ing  by  layers  will  start,  it  is  assumed  to  start  at  the  impacted  end 
of  the  crushable  material  (in  keeping  with  experimental  observations, 
except  at  certain  impact  velocities  too  low  to  be  of  interest) . 

Thus,  if  figures  1(b)  and  1(c)  were  considered  to  represent  impact 
tests,  the  location  shown  for  the  crushed  material  requires  that  the 
material  has  been  impacted  by  the  plate  or  the  plunger.  If  the 
crushable  material  had  impacted  the  landing  surface  in  figure  1(c), 
with  the  plunger  resting  on  top  of  the  material,  there  would  have  to 
be  some  crushing  at  the  landing  surface  in  conjunction  with  the 
plunger  penetration  (or  "payload  penetration"). 

9.  It  is  assumed  that  there  is  no  section  of  crushable  material  weak 
enough  to  permit  penetration  by  any  crushable  material  (i.e.,  pene¬ 
tration  by  anything  but  the  payload).  The  validity  of  this  assump¬ 
tion  is  investigated  in  a  later  section. 

10.  The  density  of  the  crushable  material  is  assumed  uniform  in  the 
uncrushed  condition. 

11.  The  compacting  strain  of  the  crushable  material  is  assumed  uniform. 

12.  Pure  vertical  translation  is  assumed. 

13.  The  analysis  neglects  all  ringing  and  focusing  of  stress  waves. 

14.  A  weightless  exterior  cover  for  the  crushable  material  is  assumed 
that  is  strong  enough  to  prevent  shattering  of  the  crushable  mate¬ 
rial.  (For  comparison  with  the  experimental  results  of  Cundall  and 
of  reference  10,  however,  the  mass  of  the  cover  employed  is  assumed 
uniformly  dispersed  throughout  the  crushable  material.) 

15.  When  there  is  payload  penetration,  that  is,  relative  motion  between 
the  payload  and  the  crushable  material,  the  two  are  assumed  perfectly 
unbonded.  This  rules  out  the  interesting  design  possibility  of  pene¬ 
tration  despite  bonding  and  means  there  is  no  need  to  include  the 
effect  of  tensile  stresses  over  the  upper  surface  of  the  payload. 

16.  It  is  assumed  that  the  landing  surface  is  perfectly  flat  and 
perfectly  rigid. 

17.  Perfect  rigidity  is  assumed  for  the  payload. 
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Summary  of  Analytical  Development 


The  general  vertically  symmetrical  landing  geometry  without  tipover  or 
horizontal  velocity  is  illustrated  in  figures  2  and  3,  which  are  used  for 
developing  the  governing  equations  and  for  the  definition  of  terms.  (Note  that 
with  gravity  terms  being  small,  "vertical'^  can  be  any  direction  that  is  both 
parallel  to  the  resultant  impact  velocity  and  normal  to  the  landing  surface.) 

In  figure  2  the  major  limiting  assumption  is  the  absence  of  shear  deformation 
(i.e.,  there  is  no  compacted  material  dragged  outboard  of  the  payload  and  no 
compacted  material  lifted  off  the  landing  surface) . 

Evaluation  of  the  stress  integrals  and  variable  masses  in  the  governing 
equations  is  greatly  facilitated  by  introducing  the  compressive  compacting 
strain  e  and  assuming  that  it  is  uniform  throughout  the  crushable  material, 
that  is. 


z  =  =  constant  (1) 

Figure  3  shows  the  general  vertically  symmetrical  landing  geometry  for  zero 
shear  deformation  and  the  assumption  of  equation  (1) .  The  latter  assumption 
is  illustrated  by  the  fact  that  hpi  is  constant  and  mci  is  a  foreshort¬ 
ened  image  of 

Figure  4  specializes  the  geometry  to  concentric  spheres  (although  most 
of  the  resulting  simplifications  would  be  realized  as  well  by  concentric 
spherical  segments).^  It  is  further  assumed  that  the  crushable  material  has  a 
specific  direction  for  maximum  normal  crushing  stress,  that  the  material  has 
been  segmented  and  oriented  to  make  that  direction  radial,  and  that  stressed 
areas  with  normals  differing  by  an  angle  a  from  that  radial  direction  feel 
normal  stresses  a  determined  by  the  law 

a  =  Cq  cos  a  ,  a  <  90°  (2) 

where  the  restriction  a  <  90°  is  required  since  the  load  has  to  be  trans¬ 
ferred  from  the  lower  to  the  upper  hemisphere  of  crushable  material,  and 
where  the  restriction  is  automatically  met  in  all  calculations  for  nonzero 
payload  radius.  (See  ref.  5  for  alternate  anisotropic  relationships.) 

Governing  equations  corresponding  to  the  summary  just  given  are 
developed  in  appendix  A.  They  are  specialized  for  various  combinations  of 
the  following  independent  assumptions: 

1.  Neglect  of  variable  mass  effect.  In  this  assumption  the  accumulation 
of  crushed  casing  material  on  the  payload  and/or  the  landing  surface 
is  neglected.  The  assumption  is  accomplished  by  setting  =  0. 

^The  spheres  were  selected  because  of  their  ability  to  absorb  impact 
energy  in  any  direction  (i.e.,  to  handle  unoriented  impacts).  Such  impacts 
may  occur  due  to  aerodynamical ly  unstable  landing  configurations,  a  strong 
laterial  wind  with  a  steep  landing  surface,  or  terminal  guidance  failure  or 
absence . 
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2.  Neglect  of  built-up  material  effect.  In  this  assumption  stresses  are 
evaluated  at  the  surface  of  infinitely  thin  sheets  of  crushed  mate¬ 
rial  rather  than  built-up  volumes.  The  assumption  is  accomplished 
by  setting  c]^  =  1 . 

3.  Neglect  of  variable  resistance  to  payload  penetration.  In  this 
assumption  the  stress  and  force  on  the  payload  retain  their  initial 
penetration  value  for  the  entire  penetration  stroke.  The  assumption 
is  accomplished  by  setting  Fpg  (e)  =  1. 


Appendix  A  is  presented  because  a  relatively  complete  theoretical 
development,  including  variable  mass,  built-up  material,  and  payload  penetra¬ 
tion,  is  needed  in  the  literature.  The  development  is  relegated  to  an  appen¬ 
dix  because  the  details  are  not  needed  to  understand  the  rest  of  the  report. 
If  the  reader  wishes  to  locate  a  specific  result  or  derivation,  he  can  refer 
to  the  various  subdivisions  of  appendix  A  listed  in  the  Table  of  Contents. 


DESIGN  PROCEDURES  FOR  SPHERICAL  GEOMETRY 


General  Design  Conditions 

The  zero-velocity  termination  conditions  defining  the  end  of  impact  are 
equations  (A46)  and  CA47) .  For  a  minimum  weight  design,  the  termination  con¬ 
ditions  should  occur  when  the  payload  or  the  compacted  material  built  up  on 
it  (in  the  case  of  payload  penetration)  touches  the  compacted  material  built 
up  on  the  landing  surface,  that  is,  when  the  sphere  of  radius  Rp  touches 
the  region  mci  in  figure  4(b)  or  when  the  two  compacted  regions  touch  in 
figure  4(c).  The  size  of  the  compacted  regions  is  based  on  the  true  material 
compacting  strain  The  quantity  cm  is  identical  to  the  e]c  of  equa¬ 

tion  (1)  except  for  the  optional  use  of  fictitious  values  for  e]c  in  the 
equations  of  motion  of  appendix  A.  If  a  margin  of  safety  is  desired,  larger 
compacted  regions  can  be  envisioned  on  the  basis  of  a  fictitious  design 
compacting  strain  called  ej,  where 


^d  ^m 


(3) 


If  L  is  defined  as  the  distance  between  the  fictitious  compacted 
regions  at  the  end  of  the  impact  stroke,  or  between  one  such  region  and  the 
payload,  then  the  design  condition  for  contact  is 


where  Rp  is  the  payload  radius  and  R  the  overall  radius.  With  e]<; 
replaced  by  e^,  it  can  be  deduced  from  equations  (A12)  and  figure  4  that 


(4) 


L  =  R  -  Rp 


^Pmax 

^d 


(5) 
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where 


is  the  maximum  absolute  payload  displacement  during  impact. 


^Pmax 

If  ^msix  is  the  corresponding  displacement  of  the  uncrushed  crushable  mate¬ 
rial,  then  qp  =  qmax  i^  absence  of  payload  penetration.  When 
■^max 

L  >  0,  the  process  is  physically  realizable,  although  not  a  minimum  weight 
design  for  the  present  spherical  geometry.  When  L  <  0,  the  process  is  not 
physically  realizable;  but  cases  for  L  <  0  may  be  recorded  in  the  process  of 
seeking  L  =  0. 


A  second  design  condition,  having  a  less  obvious  relation  to  minimum 
weight  design,  is  concerned  with  the  maximum  design  deceleration  nj^^jgg, 
where  g^  is  the  acceleration  due  to  gravity  on  earth  and  n^^j  the  maximum 
design  g  loading.  If  actual  payload  maximum  deceleration, 

the  design  condition  is 


^Pmax  < 
^md 


(6) 


For  certain  types  of  energy  absorbing  material  or  structure,  the  minimum 
weight  design  calls  for 

^Pmax  _  ^ 

^md 


(7) 


For  other  types,  any  value  of  ^Pjj^a^x^^nid  satisfying  equation  (6)  may  yield  a 

minimum  weight  design,  with  the  critical  parameter  being  a  property  of  the 
crushable  casing  (such  as  Oq)  . 


The  quantity  equations  (6)  and  (7)  is  the  maximum  of  the 

quantity  in  equations  (A32)  and  (A41) .  (See  also  eqs .  (A48) ,  (A51), 

(A55) ,  and  (A59)  for  various  specializations.)  If  variable  mass  and  built-up 
material  are  neglected  in  the  analysis  (by  setting  p^k  =  0  and  =  I, 
respectively)  ,  then  np  occurs  at  q  =  q^  =  R/2  in  the  absence  of  payload 

•^max  IT 

penetration,  providing  q  or  qp  becomes  that  large. 

A  minimum  weight  design  is  sought  for  a  given  payload  mass,  mpo, 
payload  radius,  Rp,  impact  velocity,  Uq,  design  compacting  strain,  cd?  and  max¬ 
imum  permissible  g  loading,  with  the  material  choice  being  arbitrarily 

limited  to  two  classes  of  crushable  material.  Where  equation  (7)  gives  the 
minimum  weight  design,  the  quantities  to  be  determined  are  the  overall  radius, 
R,  the  maximum  normal  stress,  Qq,  and  the  density,  of  the  crushable  mate¬ 

rial  (from  which  the  weight  can  be  calculated) .  If  the  optimization  permits 
equation  (6),  then  a  property  such  as  Oq  must  be  specified;  and  R  and  p^m 
remain  to  be  determined  (with  ^  by-product).  In  either  case,  the 

determination  requires  that  a  relation  between  Oq  and  p^^^  be  known  for  the 
crushable  material,  thereby  effectively  reducing  the  number  of  unknowns  by 
one. 
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The  required  relation  between  Oq  and  can  be  determined 

experimentally  for  a  variety  of  materials  and  stated  directly  as  in  refer¬ 
ence  2  (p.  259).  It  is  common,  however,  to  introduce  a  parameter  for  which 
values  are  widely  known,  namely,  the  specific  energy  absorption  (SEA)  which 
is  defined  as 


SEA  =  En,  I  ~ — )  (S' 

yPcmge  J 

and  which  should  generally  be  as  large  as  possible.  The  relation  between 
Oq  and  Pcm  can  be  expressed  by  giving  SEA  in  terms  of  Oq.  With  o  =  Cq 
and  e  =  e^,  the  product  OqC^  is  the  area  enclosed  by  the  dashed  lines  in 
figure  1(a),  which  is  the  energy  absorbed  per  unit  volume.  Dividing  this 
product  by  PcmSe  gives  equation  (8)  and  shows  SEA  to  be  the  energy 
absorbed  per  unit  mass . 

Figure  5  shows  the  variation  of  SEA  with  Oq  for  a  balsa- like  class 
of  material  and  a  honeycomb- like  class  of  material  (neglecting  glue- joint 
weight  and  the  effects  of  high  impact  velocity  and  low  temperature) .  The 
equation  for  SEA  in  terms  of  Oq  given  for  the  balsa- like  material  in  fig¬ 
ure  5  is 

SEA  =  24,000  ft-lb/lb  , 

SEA  =  ft-lb/lb  , 

The  balsa- like  material  is  balsa  in  the  sense  that  the  second  of  equations  (9) 
is  deduced  from  the  curve  fit  of  reference  2  (p.  259)  which  closely  approxi¬ 
mates  the  data  therein  for  soliH  balsa  of  various  densities.  The  material 
is  considered  to  be  only  balsa-like,  however,  because  the  first  of  equations 
(9)  and  the  dashed  portion  of  the  balsa-like  curve  in  figure  5  are  assumed  to 
be  valid  for  a  hypothetical  cored  balsa  (for  which  cores  of  material  are 
removed  in  the  radial  direction) .  The  assumption  is  that  roughly  one-third 
of  the  material  can  be  removed  from  the  lightest  (lowest  density)  solid  balsa, 
giving  a  Oq  range  from  800  to  1,200  psi,  without  reducing  the  SEA  by 
introducing  significant  buckling,  end  effects,  and/or  Poisson's  ratio  effects; 
and  this  assumption  seems  reasonable.  It  should  be  noted  that  the  only  cor¬ 
roborative  case  where  SEA  decreases  with  increasing  in  reference  4  is 

for  dry  balsa  (0-percent  moisture),  atmospheric  pressure,  and  an  ambient  tem¬ 
perature  of  approximately  78°  F. 

The  curve  for  the  honeycomb- like  material  is  presented  in  figure  5  so 
that  results  can  be  deduced  for  lower  SEA  values  and  also  for  the  relatively 
common  case  where  SEA  increases  with  increasing  Og.  The  equation  given  in 
figure  5  is 

SEA  =  478.5  ft-lb/lb  ,  600  psi  ^  Pq  " 


800  psi  S  =  1,200  psi 


(9) 


1,200  psi  =  Oq  =  1,800  psi 
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The  material  is  called  honeycomb -like  rather  than  a  specified  honeycomb 
because  equation  (10)  is  deduced  from  a  curve  fit  in  reference  2  (p.  259) 
that  rather  loosely  fits  data  for  several  types  of  aluminum  and  fiberglass 
honeycomb  having  a  variety  of  densities.  It  should  be  noted  that  reference  2 
makes  no  mention  as  to  whether  equations  (9)  and  (iO)  incorporate  the  dynamic 
effects  of  damping  and  dynamic  buckling  (see  item  3  in  "Summary  of 
Fundamental  Assumptions  and  Limitations") . 


Design  Procedure  for  Simplified  Model  Without  Payload  Penetration 


For  the  analytical  model  labeled  "simplified"  in  the  present  report,  it 
is  assumed  that  variable  mass  can  be  neglected  (setting  p^k  =  0) >  that  built- 
up  material  can  also  be  neglected  (c]^  =  1)  ,  that  gravity  forces  are  insignifi¬ 
cant  (g  =  0) ,  and  that  the  resistance  to  payload  penetration  is  constant 
(Fpo(e)  =  1).  Under  the  first  three  of  these  assumptions  (with  the  last 
required  only  for  penetration),  equations  (A55)  and  (A62)  apply  (in  a  simpli¬ 
fied  form  with  g  =  0)  and  the  division  of  the  latter  by  the  former  gives 
(with  definitions  from  eq.  (A54)) 


2 

V 

•'  max 


np/nmd  ^pggRp 


(R/Rp)  -.[(2/3)y^^^] 
y[(R/Rp)  -  y] 


(113 


Equation  (A55)  is  a  parabola  with  at  y  =  (R/Rp)/2.  Hence,  for  the 

present  restrictive  case, 


'max 

is  the  acceleration  at  y 


max 


as  long  as 


Xmax  ^  (R/Rp)/2.  If  ^  (R/Rp)/2,  ^p^^^g^  is  the  acceleration  at 

y  =  (R/Rp)/2.  Thus  equation  (11)  can  be  expressed  in  terms  of  np  as 


’^PmaxSeRp 


W3 

1  “  (^p/f^)ymax 


max 


1  R 


^  <^12) 
Rp 


’^Pmax^e’^p  ^  ^  '"^x/  ’ 


y 


max  2  Ri 


(13) 


If  equations  (12)  and  (13)  are  multiplied  through  by  Rp/R  and  if  the  first 

of  equations  (A54)  is  used  (giving  y  =  z  (R/R„))  the  result  is 

max  max  Jr 


z  [1  -  (2/3) z  ] 
max'-  '■  ^  max-* 


^Pmax^e*^ 


1  -  z 


max 


max 


(14) 
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(15) 


U. 


2 

3 


max 


1 

2 


Equations  (14)  and  (15),  which  are  applicable  only  in  the  absence  of 
payload  penetration,  are  plotted  in  figure  6  as  versus  (np^^^ge)P* 


Figure  6  can  be  used  to  determine 
and  payload  radius 

interfere  with  q 


max 


Zmax  ""  %ax/^  basis  of  Uq,  R, 

Rp  simply  has  to  be  small  enough  not  to 


An  interesting  feature  in  figure  6  is  that  ^ma.x  determined  with  no 
knowledge  of  Oq,  iHqq,  or  mpo-  Once  Zjj^g^x  known,  however,  m^o  can  be 
found  for  a  given  mp^  and  according  to 

2110  R^ 

O  2 

m  =  - 

CO  2  ^ 

^o 


1  “  —  z 
^  3  ^max 


-  m. 


po 


(16) 


as  derived  from  equations  (A62)  and  (A39)  with  g  =  0.  It  is  still  not 
necessary  to  know  Rp.  For  a  minimum  weight  design,  however,  Rp  must  be 
known  or  determined;  and  equation  (4)  must  be  satisfied  according  to  the 
definition  in  equation  (5).  In  addition,  the  SEA  ,  as  defined  in  equation 
(8),  may  be  known  instead  of  Oq.  When  equations  (4),  (5),  and  (8),  with 
Zmax  =  =  ^Pmax^^  "’co  "  (4/3)Trp^j^(R3  -  Rp^)  ,  are  introduced  into 

equation  (16) ,  the  result  is 


g  (SEA)  /^A/  ni  \  1  -  (z  /e,)  +  [  (z 

/  d  w  CO  \  _  ^  max  d  ^  max  d^  •  ^ 

\  in/\po  ^  ^co /  z  [(1/2)  -  (z  /3)] 

o  \  f  \  ^  max'  ^ 


where  the  symbol  is  introduced  for  convenience.  Equation  (17)  can  be 

used  to  determine  m^o  ^ov  a  given  Zj^^x  (which  implies  knowledge  of  Uq  and 
ge  as  indicated  in  fig.  6)  if  ed,  SEA,  and  mpo  are  known  (with  Cjn 
actually  canceling  the  same  quantity  in  SEA).  In  figure  7,  equation  (17)  is 
plotted  for  ej  =  0.7,  0.8,  and  0.9.  The  three  cross-plotted  values  of  Rp/R 
serve  as  a  reminder  that  figure  7  represents  a  design  for  which  contact  would 
occur  between  the  payload  and  the  compacted  material  if  that  material  had  a 
compacting  strain  of  instead  of  Cjji  (see  the  sketch,  fig.  7). 

In  the  absence  of  payload  penetration,  figures  6  and  7  are  sufficient 
for  the  simplified  model  if  R  is  given  and  Rp  is  to  be  determined.  If 
Rp  is  given,  figures  6  and  7  remain  useful  as  a  check  and  as  a  means  of 
determining  z^ax^  but  two  different  figures,  based  on  Rp,  are  more  useful 
for  the  original  design.  The  first  of  these  is  determined  by  introducing 
equations  (4)  and  (5)  into  equations  (13)  and  (12)  as  a  substitute  for 

ymax  =  Praax/Rp  =  Ppmax^^P’ 
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u 


"Pmax^e’^p 


4e^2[i  .  (Rp/R)]2 


1  -  (2/3)e^[l  -  (Rp/R)] 

"yR 


1 

(18) 


U 


Hn  g  R 

Pmax  e  p 


e^[l 


(Rp/R)]  (  1  -  (2/5)e^[l  -  (Rp/R)] 
(Rp/R)  [  1  -  e^[l  -  (Rp/R)] 


(19) 


Equations  (18)  and  (19)  are  plotted  for  =  0.7,  0.8,  and  0.9  in  figure  8. 

Figure  9  is  the  companion  to  figure  8  and  is  found  by  using  equations 
(4)  and  (5)  to  define  a  replacement  for  ^  ‘imax/^  "  '^Pmax'^^  equa¬ 
tion  (17).  The  Rp/R  found  in  figure  8  then  determines  in  figure  9, 

and  Jma  determines  m^o  as  before.  The  sketches  in  figures  8  and  9  show 
the  same  design  configuration  as  that  in  figure  7,  leaving  figure  6  as  the 
only  one  for  which  contact  is  not  required  for  compacted  material  based 
on  Cjj. 


Design  Procedure  for  Simplified  Model  With  Payload  Penetration 

Payload  penetration  was  not  permitted  in  figures  6  through  9  (because 
of  a  hypothetical  perfect  bond).  If  penetration  is  now  permitted  (because  of 
the  total  absence  "of  bonding),  the  first  step  is  to  determine  whether  it  will 
occur.  It  will  occur  if  the  design  without  penetration  gives 


where  Zg  is  the  dimensionless  displacement  at  which  the  area  of  material 
being  crushed  becomes  large  enough  to  cause  sufficient  deceleration  for  pene¬ 
tration.  The  quantity  Zg  is  defined  by  equation  (A58)  as  modified  by  the 
change  of  variable,  yg  =  Zg (R/Rp) .  The  modified  equation,  plotted  in  figure 
10,  is 


(21) 


The  curve  is  cut  off  at  Zg  =  0.5,  which  is  the  maximum  value  at  which  pene¬ 
tration  can  occur  for  the  present  approximation  (as  seen  by  the  maximum  force 
in  eq.  (A52)).  Figure  10,  while  not  ultimately  essential  for  design  purposes, 
is  useful  as  a  check,  and  the  magnitude  of  Zg  is  significant  in  evaluating 
the  importance  of  penetration  or  potential  penetration. 
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The  first  of  the  two  main  design  figures  for  penetration  is  based  partly 
on  the  first  of  equations  (A59)  and  the  last  of  equations  (A39)  with  g  =  0 
and  np  =  =  constant.  This  yields,  as  expected  for  the  present 

assumptions , 


‘Pmax  mpoge 


Equation  (22)  is  used,  together  with  equation  (21),  equation  (A64)  for  g  =  0, 
the  last  of  equations  (A39) ,  the  first  and  second  of  equations  (A54) ,  and 
equations  (4)  and  (5)  for  ^  ^Pmax'^^^  identify  the  variable  Npu  as 


follows : 


d  Pmax  e  p 


1  -  (4/3) z 


The  purpose  in  isolating  is  to  obtain  a  single  unknown  (Rp/R)  in  terms 

of  Zs  and  with  a  factor  that  is  a  function  of  known  quantities  (Uq,  ed^ 
^Pmax^  ge>  Rp)  •  The  next  step  is  to  isolate  a  variable  for  i^co^^o  'Without 

Rp/R.  This  is  accomplished  by  dividing  equation  (21)  by  the  square  of  equa^ 
tion  (23).  The  resulting  variable,  called  is 


1  +  (mco/mpo) 


2Z3CI  -  Z3) 


[1  (l  -  (Zs/Ed){[l  '  (4/3)Zs]/2(l  -  Zs)|) 


Figure  11  is  a  plot  of  N|^u  versus  for 

The  plot  is  constructed  by  selecting  numbers  for  Zg 
calculating  the  corresponding  values  of  Nru  and  N^u 
(23)  and  (24).  A  relationship  between  R  and  m^o  i 
11  in  terms  of  the  known  quantities  Uq,  ed^  ^P^ax* 

trast  to  fig.  8  without  penetration,  where  R  can  b 
quantities  and  used  in  determining  m^^Q  in  fig.  9). 
shows  the  design  condition  implicit  in  equations  (4) 
namely,  contact  between  the  two  regions  of  compacted 
replaced  by  ej* 


(24) 

=  0.7,  0.8,  and  0.9. 
between  0  and  0.5  and 
according  to  equations 
s  established  in  figure 
ge,  Rp,  and  mpo  (in  con- 

e  determined  from  known 
The  sketch  in  figure  11 
and  (5)  for  penetration, 
material  if  were 


A  definition  and  a  volume  density  relationship  are  useful  at  this  point, 
namely: 


^po  (^co/^po) 


"'co  4  F/r  V 


(p  g  )TrR^  (p  g  )7tR^  P  ttR^ 

^^cm^e^  V  ^  cm^e"^  P  cm  P 


3  3  I  R 


-  1  (25) 


where  Wpo  and  Wco  are  the  original  weights  of  the  payload  and  crushable 
material,  respectively.  Equation  (25)  is  plotted  in  figure  12;  it  is  useful 
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with  or  without  payload  penetration  and  does  not  depend  on  any  assumptions. 
For  the  case  of  payload  penetration,  however,  and  for  the  special  assumptions 
of  the  simplified  design,  equation  (22)  applies  and  can  be  substituted  into 
equation  (25)  with  equation  (8)  to  give 


^ma 


SEA 

^Pmax^^^P 


(26) 


With  equation  (26),  iHcq/^o  determined  from  Nj^^-  without  knowing  Pcm- 

Thus  figure  12  becomes  the  companion  for  figure  11  in  an  iterative  design  pro¬ 
cedure  for  payload  penetration.  The  procedure  is  simply  to  select  an  initial 
value  for  mco/n^po^  calculate  Nmu  in  terms  of  mco/^ipo  nnd  known  quantities, 
determine  Nru  from  figure  11,  evaluate  Rp/R  in  terms  of  N|^  and  known 
quantities,  determine  from  figure  12,  calculate  a  second  value  of 

mco/nipo  in  terms  of  N^a  and  known  quantities,  and  repeat  the  process  until 
two  values  of  ngree  to  the  accuracy  permitted  by  the  figures . 

The  iterative  procedure  just  described  could  have  been  avoided,  of 
course,  by  combining  equations  (21)  through  (26)  into  a  polynomial  for  z^- 
This  polynomial  could  be  solved  for  selected  parametric  values  of 
UQ^/2ejjnp^^^g^Rp  and  SE A/np^^^e^^Rp ,  together  with  values  of  such  as  those 

selected  for  the  iterative  charts,  and  Rp/R  and  mco/^^^po  could  be  determined 
accordingly.  Such  a  procedure  was  avoided,  however,  because  of  the  strong 
likelihood  that  a  reasonably  limited  group  of  parametric  values  would  not  have 
sufficiently  broad  applicability. 


Specialization  of  Design  Procedures  for  Simplified  Model 
to  Materials  of  Figure  5 


Regardless  of  whether  penetration  is  absent  or  present,  it  is  apparent 
that  all  of  the  design  figures  (6  through  12)  are  essentially  independent  of 
the  material  or  structure  selected  for  energy  absorption.  To  calculate  m^o 
from  figures  7,  9,  and  12,  however,  an  SEA  has  to  be  selected.  With  mco 
and  R  determined  by  the  figures  and  Rp  known  in  advance,  can  be  calcu¬ 

lated  (with  the  aid  of  fig.  12  if  desired);  and  Oq  can  then  be  found  for  the 
selected  SEA  according  to  equation  (8) .  All  this  implies  the  assumption 
that  a  material  having  the  calculated  properties  is  available  or  that  a 
corresponding  structure  can  be  constructed. 


For  the  comparison  purposes  of  this  paper,  however,  the  materials  are 
restricted  to  a  choice  between  the  two  classes  described  in  figure  5.  With 
figures  6  through  12  having  been  constructed  independently  of  figure  5  (in 
the  interest  of  generality) ,  the  use  of  figure  5  imposes  a  trial-and-error  (or 
transcendental)  solution  for  those  designs  in  which  payload  penetration  is 
absent.  The  trail-and-error  solution  is  aided  by  incorporating  equations  (8) 
and  (25)  into  equation  (17)  to  yield 


SEA  = 


^m*^ma 


(Se^d/^o^)  "  C'^ma^po^^o^ma^^p^) 


(27) 
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On  the  basis  of  given  values  of  em,  ge>  ^d,  ^po?  ^Pmax" 

values  of  Jma  ^ma  determined  according  to  figures  8,  9,  and  12,  the 
quantities  emdina^  ge^d/^o^^  and  dma^po/^ma'^Rp^  in  equation  (27)  can  be 
calculated  in  advance.  Equation  (27)  and  figure  5  then  become  the  basis  for 
the  numerical  solution. 


A  sample  calculation  for  the  trial- and- error  solution  just  described  is 
given  in  appendix  B.  It  should  be  noted  that  the  trail  and  error  of  appendix 
B  would  have  been  eliminated  if  the  SEA  had  been  constant  -  despite  varia¬ 
tions  in  Go,  as  for  the  balsa-like  curve  in  figure  5  when  Oq  <  1,200  psi. 
There  is,  however,  another  type  of  trial  and  error  that  occurs  when  a  spe¬ 
cific  value  of  On,  rather  than  n^  ,  is  sought.  In  this  case,  calcula- 

Fmax^  ^ 


tions  like  those  in  appendix  B  (but  without  the  Oq 
performed  for  successive  selections  of  until 
desired  value. 


-  SEA  trials)  must  be 
Oq  approaches  the 


If  (in  contrast  to  appendix  B)  payload  penetration  is  present,  the  use  of 
the  specific  materials  in  figure  5  does  not  impose  a  trail-and-error  solu¬ 
tion.  The  reason  is  that  the  unknown  m^o  ^oes  not  appear  in  equation  (22) . 
In  addition,  the  simplicity  of  equation  (22)  means  that  either  Oq  or 
can  be  selected  without  requiring  trial  and  error. 

A  sample  calculation  for  the  iterative  penetration  procedure  described 
earlier  is  given  in  appendix  C  for  the  balsa- like  material  of  figure  5. 

Only  three  iterations  are  needed  in  appendix  C  because  of  a  fortunate  initial 
guess  of  mco/nipo-  Th^  initial  guess  was  based  for  all  calculations  on  other 
penetration  cases  or  prior  examples  without  penetration,  and  the  worst  guess 
required  five  iterations.  When  a  penetration  design  is  not  feasible,  either 
the  Oq  value  of  equation  (22)  will  be  beyond  the  range  of  figure  5,  or  the 
iterations  will  move  off  the  curves  of  figures  11  and  12, 


Computer  Procedures 

A  set  of  three  computer  procedures  has  been  programmed  to  evaluate  the 
governing  equations  of  the  impact  problem  for  spheres  (eqs .  (A32) ,  (A37) ,  and 
(A41))  in  accordance  with  the  appropriate  termination  conditions  (eqs.  (A46) 
and  (A47))  and  the  appropriate  design  conditions  (eqs.  (4),  (6),  and  (7)). 
These  procedures  are  described  in  appendix  D.  They  can  be  used  to  check  the 
simplified  model  designs  or  to  initiate  simplified  designs  and  other  designs 
based  on  more  complicated  models. 

Specifically,  of  the  three  computer  procedures,  two  are  designs  in  that 
they  automatically  iterate  initial  guesses  to  determine  required  crushable 
casing  parameters.  One  of  the  two  design  procedures  varies  the  overall  radius 
R  and  the  material  maximum  crushing  stress  to  achieve  a  desired  accel¬ 

eration,  as  indicated  by  equation  (6)  or  equation  (7),  and  a  desired  ratio  of 
stroke  to  available  stroke,  as  suggested  by  equation  (4).  This  program  is 
applicable  only  in  the  absence  of  payload  penetration.  The  other  design 
varies  only  R  for  a  selected  Oq  (fixed  material  once  SEA  or  a  plot  of 
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SEA  versus  Oq  is  selected)  and  achieves  only  the  desired  stroke  ratio.  It 
applies  both  with  and  without  penetration. 


Penetration  requires  only  a  search  for  R  since  Oq  is  essentially 
determined  in  that  case  regardless  of  R  or  the  stroke.  The  reason  becomes 
apparent  when  the  first  of  equations  (A41)  is  modified  at  the  start  of  pene¬ 
tration  (e  =  de/dx  =  0)  by  equations  (A45) ,  (A42) ,  and  the  last  of  equations 
(A39)  to  give  the  expected  result 


np 


max 


TlRp^Oo 

mpoSe 


(28) 


where  the  approximate  equality  sign  is  used  because  n„ 

Pmax 

from  the  value  at  the  start  of  penetration  for  the  detailed 

one  part  in  a  thousand  for  the  present  examples) .  Equation 

Oq  for  a  given  n„  regardless  of  R  or  the  stroke. 
Pmax 


differs  slightly 

model  (by  at  most 
(28)  determines 


The  third  computer  procedure  can  also  be  used  with  or  without  penetration 
but  is  not  programmed  to  iterate  and  produce  a  design.  Hence  it  can  be  used 
only  to  check  a  given  configuration. 


It  should  be  noted  that  all  three  programs  permit  SEA  to  be  selected 
arbitrarily  or  calculated  (after  Oq  is  selected)  according  to  equations  (9) 
and  (10)  for  the  materials  considered  herein. 


Design  Procedure  for  Detailed  Model  With  and 
Without  Payload  Penetration 

The  detailed  model  is  the  second  of  the  two  approximate  analytical 
models.  Where  the  simplified  model  had  Pck  =  0,  =  1,  g  =  0,  and 

Fpo(e)  =  1,  the  detailed  model  has  =  Pcm>  ^k  =  ^m^  8  =  gL>  and  from 

equation  CA45) , 

+  l)Sp  dSp 


p  +  e^  +  1 

The  Pq]^  =  equation  means  that  the  detailed  model  incorporates  variable 

mass  (accumulating  on  the  payload  and/or  the  landing  surface)  according  to 
the  crushable  material  density,  and  implies  a  finite  volume  of 

compacted  material  instead  of  an  infinii:ely  thin  sheet  to  determine  the  sur~ 
face  for  stress  evaluation.  The  equation  g  =  gL  simply  incorporates  a 
negligibly  small  gravity  term  for  completeness,  and  the  integral  for  Fpo(e) 
permits  a  calculated  deviation  from  a  constant  resistance  to  penetration, 

l^ile  the  computer  procedures  described  in  the  previous  subsection 
constitute  an  alternative  to  figures  5  through  12  for  the  simplified  model, 
they  constitute  the  only  design  procedure  presented  herein  for  the  detailed 
model.  Of  course,  figures  5  through  12  remain  useful  as  starting  points  for 
the  computer  iterations . 


F  fe) 
po 


(es. 


/2^ 
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DESIGN  RESULTS  AND  DISCUSSION  FOR  SIMPLIFIED  AND  DETAILED 
MODELS  HAVING  SPHERICAL  GEOMETRY 


Description  of  Design  Examples  and  Most  General  Results 

Thirteen  examples  have  been  calculated  by  the  simplified  and  detailed 
analytical  models  described  in  the  previous  section.  The  impact  velocity 
Uq  is  taken  to  be  300  ft/sec  for  all  examples,  the  material  compacting 
strain  sm  to  be  0.8,  the  design  compacting  strain  ed  to  be  0.7,  and  the 
payload  maximum  g  loading  to  be  2000  or  less.  The  examples  vary 

in  the  given  payload  weight  Wpo.  the  given  payload  radius  Rp,  the  presence 

or  absence  of  a  hypothetical  perfect  bonding,  and/or  the  given  material. 

All  examples  are  based  on  a  choice  between  two  classes  of  crushable 

material  -  the  balsa- like  and  honeycomb- like  materials  defined  in  equations 

C9)  and  flOl  and  shown  in  figure  5.  The  desired  value  of  n^  (for  which 
^  ^max 

the  label  n^^s  used  where  needed)  is  given  for  some  examples  and  the 

maximum  crushing  stress  Oq  for  others. 

The  results  for  all  examples  are  presented  in  table  1  in  terms  of  the 
presence  or  absence  of  penetration,  the  value  of  ^p^^^  (^^  given)  or  of 

Oq  (if  not  given),  the  crushable  material  density  PcmSe^  specific  energy 

absorption  SEA,  the  stroke-to-potential-start-of-penetration  ratio  9pj^^^/9s> 

the  overall  radius  R,  the  crushable  material  weight  the  total  weight 

Wco  +  dimensionless  unused  stroke  L/R.  When  L/R  is  negative, 

the  payload  would  go  too  far  (slightly,  for  the  present  examples)  and  cause 
excessive  accelerations  if  it  were  not  for  the  stroke  margin  of  safety  given 
by  the  use  (in  defining  L/R)  of  the  fictitious  design  compacting  strain  ed 
rather  than  the  material  value  Cm* 

The  resulting  quantities  just  listed  are  presented  for  the  simplified 
and  detailed  models  in  table  1;  the  corresponding  ratios  of  detailed  to  sim¬ 
plified  results  are  presented  when  the  results  are  numerical.  It  is  apparent 
from  the  W^o  ratios  that  the  simplified  model  has  the  lower 

An  important  ground  rule  for  the  ratios  of  table  1  is  that  either  the 

On  ratio  or  the  n^  ratio  is  required  to  be  essentially  unity.  The 

^max 

choice  is  made,  after  the  simplified  model  has  been  calculated,  in  favor  of 
the  lightest  resulting  detailed  model.  This  gives  the  best  weight  ratio, 
that  is,  the  ratio  closest  to  unity. 

In  view  of  this  choice,  it  is  not  surprising  that  the  worst  Oq  or  n^ 

^max 

ratios  (of  detailed  to  simplified  results)  are  farther  from  unity  than  the 
worst  weight  ratios.  In  fact,  the  worst  ratio  of  all  is  1.2917  for  Oq  in 
case  1  (except  for  L/R,  where  numbers  approach  zero,  fortunately).  This 
ratio,  although  not  excessively  different  from  unity,  represents  a  large 
enough  difference  in  required  material  (within  the  honeycomb -like  category  of 
case  1)  that  a  designer  would  wish  to  know  whether  the  simplified  model  is 
numerically  more  realistic  than  the  detailed  model  or  vice  versa. 
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In  a  footnote  of  table  1  it  is  noted  that  the  simplified  results  are 
*^chart”  results  (derived  from  figs.  5  through  12)  except  for  the  L/R  column 
and  cases  7  and  12,  which  are  automatic  computer  results.  The  two  cases 
(both  without  penetration)  were  added  after  the  computer  design  programs 
became  available  and  were  calculated  by  those  programs  for  convenience. 

The  simplified  results  determined  by  the  chart  method  have  been  checked 
by  the  computer  checking  program  (as  opposed  to  a  design  program) .  The 

derived  computer  quantities  n^  and  /qc  had  a  maximum  error  of  1 

^max  -^max 

percent  relative  to  the  chart  method,  and  the  derived  computer  quantity  L/R 
was  close  to  zero  (ranging  from  -0.001781  to  0.001163  with  an  average  of 
-0.000194)  as  compared  to  a  zero  chart  value.  Slightly  larger  differences 
between  the  chart  method  and  the  computer  checking  method  occurred  for 
and  +  WpQ,  but  these  differences  resulted  purely  from  errors  in  calcu¬ 

lating  crushable  material  volumes  in  the  chart  method.  The  chart  and  compu¬ 
ter  results  combined  to  form  smooth  plots  for  the  simplified  model,  as  seen 
in  figures  13  and  14. 

All  numerical  results  from  table  1  are  plotted  in  figures  13  and  14 

except  for  +  Wpo  (which  is  considered  less  important  than  ,  L/R 

(which  is  often  sporadically  variable  in  sign),  cases  3,  5,  and  the  detailed 

model  for  case  2  (which  will  be  discussed  later) .  Note  that  the  curves  for 

the  simplified  and  detailed  models  are  roughly  parallel  (and  only  slightly 

curved)  over  their  mutual  abscissa  range.  Hence  trends  are  the  same  for  the 

two  models.  Even  the  crossover  in  figure  14(d)  for  q  /qc  without  pay- 

Pmax 

load  penetration  is  almost  parallel  and  thus  maintains  the  trend  with  only  a 
slight  comparison  reversal  between  models. 


Effect  of  Payload  Radius  (and  Payload  Packaging  Density) 
for  Balsa- like  and  Honeycomb -like  Materials 
Without  Payload  Penetration 

Crushable  casing  properties  and  performance  are  presented  in  figure  13 

(with  and  without  penetration)  as  functions  of  payload  radius  Rp  (with  four 

numbers  attached  for  payload  packaging  density  Ppge)  honeycomb- like 

material,  a  payload  weight  of  100  lb,  and  an  n.^  value  of  2000.  The  cor- 

Pmax 

responding  cases  in  table  1  are  1,  4,  6,  and  7  without  penetration  and  case  2 
(simplified  model  only)  with  penetration. 

In  the  absence  of  payload  penetration,  figures  13(a),  13(b),  and  13(c) 
show  that  an  increasing  Rp  (decreasing  Ppge)  causes  SEA,  R,  Oq,  and 

PcmSe  increase  (with  held  at  the  maximum  permissible  value  of 

2000,  which  minimizes  according  to  preliminary  calculations).  The  poten¬ 
tial  penetration  ratio  exists  only  for  the  simplified  model  at 

Rp  =  0.6  ft  (among  the  plotted  points),  indicating  penetration  to  be  impossible 
for  the  other  cases.  Between  the  two  models,  the  simplified  model  has  the 
lower  (slightly),  the  lower  SEA,  the  higher  R,  the  lower  Oq,  and  the 

lower  PcmSe*  The  lower  Pcmge  the  only  apparent  reason  for  the  lower 
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W^Q.  It  results  from  the  lower  which  is  made  possible  by  the  fact  that 
the  simplified  model  sees  a  higher  average  stress  for  a  given  Oq  than  does 
the  detailed  model. 


Somewhat  similar  curves  are  shown  in  figure  14  for  balsa-like  material 
and  a  payload  weight  of  450  lb  (the  cases  from  table  1  being  8,  10,  12,  and  13 
without  penetration,  and  9  and  11  with  penetration).  It  should  be  noted  that 
the  450- lb  payload,  in  contrast  to  the  lighter  payload,  permits  the  use  of  the 
efficient  balsa-like  material  without  causing  excessive  decelerations  in  the 
absence  of  penetration.  When  penetration  is  absent,  results  are  presented  for 
the  lowest  Oq  considered,  namely,  800  psi.  According  to  preliminary  calcula¬ 
tions,  this  Oq  value  gives  minimum  weight  results  for  the  balsa- like  material 

(cored  when  is  less  than  1200  psi)  ,  as  opposed  to  n^^  =  2000  for  the 

^max 

honeycomb- like  material.  The  stress  of  800  psi  determines  SEA  as  24,000 

ft-lb/lb  and  PcmSe  3,84  Ib/ft^,  Figures  14(a),  14(b),  14(c),  and  14(d) 

then  show  that  an  increasing  Rp  (decreasing  Ppge)  causes  R,  and 

nr)  to  increase  but  causes  /q^  to  decrease  (decreasing  the  likeli- 

Fmax  Pmax  ^ 

hood  of  penetration,  as  expected  for  increasing  Rp) .  Note  that  the 

9p  /9s  0.730  at  R^  =  1.6  ft  is  the  only  case  where  penetration  is 


max  P 

impossible  because  qrv  /q^.  is  less  than  one.  A  comparison  of  the  simpli- 

^  max  ^ 

fied  and  detailed  models  indicates  that  the  former  has  the  lower  the 


The 


curves  have  a  shallow  cross - 


lower  R,  and  the  higher  np 
over  (in  the  absence  of  payload  penetration) ,  and  q  /q^  does  not  exist 


max 


for  the  detailed  model  at  R, 


1.6  ft 


Pmax 


This  absence  of  q^^  /q^  for  the  detailed  model  when  it  exists  for  the 

Pmax 


simplified  model  has  been  a  recurring  theme  without  penetration,  as  seen  in 
table  1,  and  indicates  that  the  detailed  model  is  the  least  susceptible  to 
penetration.  The  reason  for  this  is  the  relatively  low  maximum  g  loadings. 


n 


Pmax 


for  the  detailed  model  in  figure  14(c). 


Effect  of  Payload  Penetration 

The  data  for  payload  penetration  in  figures  13  and  14  are  limited  by  the 
fact  that  penetration  often  does  not  occur  even  when  the  payload  is  unbonded 
(see  the  qp  ^_/q^  column  in  table  1  for  case  2  with  the  detailed  model  and 


increases  SEA  and  Oq, 

/q^  for  the  simplified 
max  ^ 


max 

for  cases  4,  6,  7,  and  13  with  both  models).  It  is  seen  in  figure  13  for 
honeycomb- like  material,  that  penetration  reduces 
decreases  R,  increases  PcmSe^  decreases  qp 

model  at  Rp  =  0.6  ft.  The  impossibility  of  penetration  for  the  plotted 
abscissas  higher  than  Rp  =  0.6  ft  means  that  the  curve  for  the  higher 

Rp  values  could  have  been  combined  with  the  penetration  point  at  Rp  =  0.6  ft 
to  form  a  curve  for  an  unbonded  payload,  with  an  obviously  beneficial 

effect. 


The  W^Q  benefit  due  to  penetration  at 
increase  in  SEA  and  Oq  and  the  corresponding  decrease  in 


Rp  =  0.6  ft  results  from  the 
R.  These 
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quantities  can  change  because  the  honeycomb- like  material  is  allowed  to  change 

within  its  class  to  maintain  =  2000  for  penetration  (which  gives  the 

^max 

lowest  Wcor  according  to  preliminary  calculations,  as  it  did  without 
penetration) . 


In  figures  14(a)  and  14(d)  for  balsa-like  material,  it  is  apparent  that 
payload  penetration  increases  and  decreases  both  models 

at  Rp  =  1,0  ft  and  Rp  =  1.2  ft  when  SEA,  Qq,  and  PcmSe  held  constant  at 

24,000  ft-lb/lb,  1,200  psi,  and  5.76  Ib/ft^,  respectively  (where  Oq  =  1,200 
psi  defines  the  elbow  of  the  SEA  curve  in  fig.  5  and  gives  the  minimum 
according  to  preliminary  calculations) .  The  overall  radius  R  is  shown  in 


figure  14(b)  to  be  increased  by  penetration  at  Rp  =  1.0  ft  and  decreased  at 
Rp  =  1.2  ft;  the  apparent  contradiction  between  the  R  and  effects  of 

penetration  is  resolved  by  recalling  that  the  lowest  Wco  without  penetration 
was  for  a  lower  PcmSe>  namely  3.84  Ib/ft^,  than  for  penetration.  In  figure 


14(c)  the  effect  of  penetration  on 


is  seen  to  be  the  same  for  the  two 


models  at  Rp  =  1.0  ft  but  different  at  Rp  =  1.2  ft. 


The  W^Q  effect  of  penetration  indicates  that  a  perfect  payload  bonding 
would  be  desirable  for  the  cored  balsa-like  material  at  Rp  =  1.0  ft  and 
Rp  =  1,2  ft.  In  the  event  that  such  a  bonding  is  not  feasible  or  trustworthy, 
however,  the  case  of  an  unbonded  payload  must  be  considered.  Hence  undefined, 
wavy- line  transitions  between  penetration  and  no  penetration  are  shown  in 
figures  14(a),  14(b),  and  14(c),  but  not  in  figure  14(d)  (due  to  lack  of 
space).  The  most  important  of  these  transitions  is  for  Wqq  in  figure  14(a). 
If  the  Wco  transition  were  specified  in  the  area  of  the  wavy  line,  it  would 
define  a  curve  for  an  unbonded  payload  since  the  wavy  line  skips  over 
Rp  =  1.4  ft,  where  bonding  is  required  to  prevent  penetration  at 
Oq  =  800  psi. 


Stresses  between  800  and  1200  psi  (the  minimum  weight  value  for 
penetration)  will  presumably  be  useful  in  the  transition.  Thus  a  variety  of 
specified  transitions  will  be  possible.  For  an  unbonded  payload,  it  would  be 
desirable  to  seek  a  minimum  in  the  transition  region  for  a  single  Oq 

at  which  penetration  is  ready  to  begin  at  the  end  of  the  stroke;  but  this  is 
beyond  the  scope  of  the  present  report.  The  reasonable  assumption,  however, 
is  that  such  a  minimum  W^o  exists  and  indicates  an  important  design  trade¬ 
off  for  cored  balsa,  a  point  where  a  further  increase  in  payload  packaging 
density  (the  abscissa  of  fig.  14(a))  is  undesirable. 


The  most  important  effects  of  payload  penetration  discussed  so  far  have 
been  a  decrease  of  for  the  honeycomb -like  material  and  an  increase  of 

Wco  the  balsa-like  material.  Both  comparisons  have  been  based  on  mate¬ 

rial  variations  for  minimum  weight  within  the  categories.  If  these  varia¬ 
tions  are  not  allowed,  penetration  can  always  be  expected  to  increase  Wcq* 

Except  for  the  roughly  7-percent  decrease  in  Wco  honeycomb¬ 

like  material  at  Rp  =  0.6  ft,  there  has  been  no  advantage  of  penetration 
reported  up  to  this  point.  If,  however,  a  design  is  permitted  to  change 
from  the  honeycomb- like  category  to  the  balsa-like  categoiy,  then  a  major 
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^co  ^po* 


penetration  advantage  occurs  in  the  present  examples  for  the  100-lb  payload. 
Specifically,  a  comparison  of  case  1  (Rp  =  0.6  ft  with  “  2000, 

honeycomb- like,  no  penetration)  and  case  3.  (Rp  =0.6  ft,  Oq  =  1200  psi, 
balsa- like,  penetration)  from  table  1  shows  a  weight  saving  due  to  penetra 
tion  by  a  factor  greater  than  3  in  and  greater  than  2  in 

In  addition,  a  comparison  of  case  4  (Rp  =  0.7  ft  with  ”  2000, 

honeycomb- like,  no  penetration)  and  case  5  (Rp  =  0.7  ft,  Oq  =  902.4  psi  at 
np^ax  of  2000,  balsa-like,  penetration)  shows  a  weight  saving  due  to 

penetration  by  a  factor  greater  than  4  in  and  almost  3  in  W^o  ^po 

Presumably,  a  still  more  impressive  comparison  would  have  occurred  at 
Rp  =  0.8  ft  except  that  a  penetration  design  could  not  be  achieved  for 


np 


max 


^  2000. 


The  honeycomb- like  material  is  used  in  the  comparisons  when  penetration 
is  absent.  The  reason  is  that  the  balsa- like  material,  having  the  higher 
SEA  over  the  same  stress  range  (fig.  5),  is  considerably  lighter  than  the 
honeycomb- like  material.  In  fact,  with  the  force  determined  by  the  stress 
range,  the  balsa-Iike  material  provides  so  little  mass  that  the  g  loading 
exceeds  the  limit  of  2000  for  the  100- lb  payload;  and  the  heavier  honeycomb¬ 
like  material  is  required  in  the  absence  of  payload  penetration.  (Note  that 
this  is  not  true  for  the  450- lb  payload.) 


When  penetration  is  permitted,  however,  for  the  balsa  design  by  removal 
of  the  hypothetical  perfect  bonding  between  the  top  of  the  payload  and  the 
crushable  material,  then  the  only  stresses  acting  on  the  payload  are  the 
crushing  stresses  at  the  bottom  of  the  payload.  This  reduces  the  force  suf¬ 
ficiently  that  the  payload  maximum  g  loading  can  be  held  to  2000  even  for 
the  balsa-like  material  with  the  100- lb  payload.  The  high  SEA  of  the  balsa¬ 
like  material  then  produces  the  large  weight  saving  due  to  penetration. 

Note  that  this  weight  saving  applies  for  the  selected  classes  of 
material;  it  would  obviously  be  decreased  if  the  honeycomb -like  material  were 
replaced  by  an  intermediate  class  just  heavy  enough  to  bring  np  down  to 

2000  without  penetration.  The  possible  availability  of  such  a  class  is 
indicated  by  the  less  efficient  balsa  reported  in  references  4  and  5.  These 
balsas  may  be  made  still  less  efficient,  in  the  sense  of  being  heavier,  by 
glue  joints  or  by  the  addition  of  weights  (a  required  weight  being  the  cover 
for  the  crushable  material)  . 


Thus  the  large  weight  saving  due  to  penetration  is  clearly  restricted 
to  the  presently  selected  materials.  It  is  obviously,  however,  a  phenomenon 
worth  considering.  The  fact  that  a  heavy  honeycomb  material  is  still  being 
considered  in  recent  design  studies  is  indicated  in  references  5,  8,  and  9; 
and  the  g  loading  issue  raised  herein  is  emphasized  in  reference  5. 


Comparison  With  Previous  Analytical  Models 

The  foregoing  results  apply  to  the  specific  analytical  models 
considered  here;  and  a  question  arises  as  to  their  validity.  They  can  be 
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considered  valid  in  one  sense  if  they  constitute  a  logical  analytical 
extension  of  a  reasonably  standard  body  of  theory  having  at  least  a  limited 
experimental  verification.  The  extended  analytical  models  can  then  be  compared 
with  future  experiments  to  establish  experience  factors  as  well  as  means  of 
improving  the  models. 

The  logical  extension  of  standard  theory  is  considered  first,  with  the 
limited  experimental  verification  to  follow.  As  pointed  out  in  the  Introduc¬ 
tion,  the  present  analytical  models  extend  a  significant  amount  of  prior  work 
(e.g.,  refs.  1  through  9)  to  include  the  effect  of  payload  penetration.  The 
summary  of  assumptions  and  limitations  given  under  ’’Outline  of  Theory"  is  the 
same  as  the  assumptions  in  references  1  through  9,  where  stated,  except  for  a 
few  variations  on  the  assumption  of  equation  (2).^ 

With  the  fundamental  assumptions  the  same,  it  is  no  surprise  that  the 
basic  equation  of  motion  (with  the  anisotropy  relation  removed)  is  the  same 
for  the  present  simplified  model  without  penetration  as  for  comparable  models 
(mass  assumed  constant,  infinitely  thin  disk  of  crushed  material,  gravity 
neglected)  in  references  1  through  9,  where  stated.  For  the  present  detailed 
model,  the  most  nearly  comparable  and  completely  described  prior  model  is 
given  in  reference  2.  The  two  models  are  identical  except  that  the  model  of 
reference  2  has  a  uniform  and  isotropic  crushing  stress  instead  of  the  par¬ 
ticular  radial  distribution  of  equation  (2) .  In  appendix  E,  it  is  shown  that 
such  a  crushing  stress  converts  the  fundamental  equation  (A19)  into  equation 
(E7),  which  agrees  exactly  with  the  fundamental  equation  (1-4)  in  appendix  A 
of  reference  2. 

The  analytical  extension  to  include  payload  penetration  is  logical  for 
the  detailed  model  in  the  sense  that  no  fundamental  assumptions  are  added  to 
the  list  given  under  "Outline  of  Theory."  For  the  simplified  model,  the 
assumption  is  added  that  the  force  resisting  penetration  is  constant,  but  it 
has  already  been  pointed  out  that  this  is  true  for  the  detailed  model  to  one 
part  in  a  thousand  or  less,  for  the  examples  considered  herein. 


Comparison  With  Previous  Experiment 

The  applicable  experimental  information  known  to  the  author  is  limited 
to  two  configurations,  one  tested  at  an  impact  velocity  of  374  ft/sec  (ref. 

10)  and  the  other  at  roughly  220  ft/sec  (footnote  1) .  The  latter  velocity 
actually  represents  an  average  of  velocities  ranging  from  215  to  225  ft/sec 
for  four  nearly  identical  tests.  At  the  higher  impact  velocity,  only  a 
single  test  is  considered  because  the  exterior  cover  for  the  model  in  that 
test  was  the  only  one  (out  of  four)  sufficiently  strong  and  resilient  to 

prevent  model  disintegration. _ 

3These  variations  in  anisotropy  can  be  important  for  various  crushable 
materials,  as  pointed  out  in  reference  5.  It  may  be  desirable  to  consider  the 
standard  theory  as  incorporating  the  variations  of  reference  5,  with  the 
assumption  of  equation  (2)  considered  as  an  example  for  evaluating  the  effect 
of  penetration.  If,  however,  it  is  desired  to  have  only  one  anisotropy  rela¬ 
tion  for  a  standard  theory,  then  the  assumption  of  equation  (2)  is  desirable 
since  it  is  simple  and  widely  used. 
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The  tests  just  described  demonstrated  the  existence  of  payload 
penetration,  depending  on  the  efficiency  of  the  payload  bonding.  They  also 
indicated  that  both  the  detailed  and  simplified  models  provide  reasonable 
(within  10  to  30  percent)  estimates  of  measured  impact  deceleration  and 
stroke.  There  are  not  enough  data  to  determine  which  model  is  quantitatively 
better. 


Alternate  Models 

The  analytical  models  considered  here  are  manifestly  only  two  out  of 
many  possibilities.  Even  with  attention  restricted  to  the  present  computer 
programs,  it  would  be  possible  to  investigate  separately  the  effects  of  the 
limitations  calling  for  Pck  ~  ^  ^k  =  rather  than  both  together  (the 
assumptions  of  g  =  0  and  Fpo(e)  =  1  being  trival) .  It  would  also  be  possi¬ 
ble  to  set  =  0  in  the  present  programs  for  integration  of  the  equation  of 

motion  but  not  for  evaluation  of  the  acceleration.  The  programs  could  be  mod¬ 
ified  to  incorporate  various  levels  of  bonding  (rather  than  just  zero  or  per¬ 
fect  bonding),  various  levels  of  shear  resistance,  and  rational  methods  of 
incorporating  the  weight  and  strength  of  the  glue  joints  and  exterior  cover. 
(Note  that  any  SEA  changes  due  to  glue  joints  or  freezing  of  the  crushable 
material  could  be  incorporated  by  changing  fig.  5  and  the  corresponding 
equations  defining  the  material.) 

Two  alternate  models  have  been  briefly  investigated;  a  full  investigation 
would  require  changes  in  the  machine  programs  or  the  design  charts.  In  one, 
called  the  ^'hemisphere"  model,  the  payload  and  the  entire  upper  hemisphere  are 
able  to  penetrate  as  a  unit,  regardless  of  bonding,  because  of  low  resistance 
to  cross-grain  crushing  in  the  equatorial  plane.  In  the  other,  called  the 
"shear-plug"  model,  the  payload  and  the  cylinder  (with  rounded  ends)  directly 
above  it  are  able  to  penetrate  as  a  unit,  again  regardless  of  bonding,  because 
of  low  shear  resistance  over  the  cylinder  walls. 

The  investigation  of  the  hemisphere  model  indicates  that  the  payload  and 
upper  hemisphere  would  have  started  to  penetrate  below  the  maximum  g  load¬ 
ings  for  almost  all  cases  in  table  1,  and  would  have  done  so  more  readily  than 
the  isolated  payload  for  the  detailed  and  simplified  models,  provided  that  the 
cross -grain  crushing  strength  is  8  percent  of  the  end  grain  value.  Hemisphere 
penetration  is  calculated  even  at  18  percent  (ref.  3  for  balsa)  for  several 
cases.  At  8  and  14  percent,  however,  it  is  also  calculated  for  the  high-speed 
(374  ft/sec)  experimental  configuration  at  a  g  loading  of  3000  to  4000  (ref. 
10);  yet  there  was  no  evidence  of  hemisphere  penetration  (i.e.,  of  exterior 
cover  wrinkling  at  the  equator).  Thus,  18  percent  may  be  the  most  nearly  cor¬ 
rect  figure,  particularly  when  the  glue  joints  are  considered.  In  any  event, 
the  hemisphere  model  is  of  interest. 

For  the  shear-plug  model,  on  the  other  hand,  there  is  no  possibility  of 
shear-plug  penetration  at  the  g  loadings  under  consideration  according  to 
calculations  based  on  a  measured  shear  strength.  This  model  would  be  of 
interest  only  for  much  higher  g  loadings  (say,  greater  than  3000)  and/or 
much  lower  shear  strengths  (less  than  145  psi) . 
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Finally,  a  question  remains  as  to  the  effect  of  payload  radius  Rp  when 
densities  (including  payload  packaging  density)  and  stresses  remain  the  same, 
that  is,  a  question  as  to  the  effect  of  scaling.  Except  for  the  trivial  grav¬ 
ity  force,  the  equations  of  motion  indicate  (for  a  nondimensionalization  dif¬ 
ferent  from  that  of  appendix  A)  that  geometric  scaling  should  apply,  with  Uq 
and  SEA  unaffected,  with  ^  1/Rp>  R  ~  Rp,  and  W^o  ~  Rp^.  These  scaling 

conditions,  then,  can  be  used  to  extend  the  present  results.  They  can  also  be 
used  to  check  numerical  results  whenever  two  cases  scale  each  other.  In  the 
present  examples,  only  cases  3  and  9  should  roughly  scale  each  other,  as 
indeed  they  do. 


CONCLUDING  REMARKS 


The  governing  equations  for  the  landing  impact  of  a  rigid  payload 
protected  by  a  crushable  casing,  including  the  possibility  of  penetration  by 
the  payload  into  the  casing,  have  been  developed  for  a  general  vertically 
symmetrical  landing.  The  general  equations  have  been  specialized  for  zero 
shear  resistance,  a  constant  compacting  strain,  and  uniform  density  of  the 
crushable  casing.  They  have  also  been  specialized  for  a  spherical  payload  and 
casing,  with  the  latter  having  its  highest  crushing  stress  in  the  radial  direc 
tion.  For  the  spherical  configuration,  two  approximate  analytical  models  have 
been  defined:  (1)  a  detailed  model  with  no  additional  assumptions  but  requir¬ 
ing  an  automatic  computer  program;  and  (2)  a  simplified  model  utilizing  either 
the  computer  program  or  design  charts  but  requiring  the  assumptions  of  infi¬ 
nitely  thin  sheets  of  crushed  material,  constant  mass  in  the  equations  of 
motion,  zero  acceleration  due  to  gravity,  and  a  constant  resistance  to 
payload  penetration.  When  specialized  to  prevent  penetration,  the  simplified 
model  (or  slight  variations  thereof)  has  been  widely  employed  in  previous 
work;  and  the  specialization  of  the  detailed  model  for  a  uniform  and  isotropic 
crushing  stress  without  penetration  has  been  shown  to  have  the  same  basic  equa 
tion  as  an  earlier  model.  Results  for  the  two  models  are  in  reasonable  agree¬ 
ment  with  two  previous  measurements,  having  impact  velocities  of  220  and 
374  ft/sec. 

Thirteen  examples  have  been  presented  for  each  of  the  two  analytical 
models.  The  examples  are  for  an  impact  velocity  of  300  ft/sec,  a  maximum 
permissible  g  loading  of  2000,  payload  weights  of  100  and  450  lb,  payload 
radii  ranging  from  0.6  to  1.6  ft,  payload  packaging  densities  ranging  from 
26.23  to  110.5  Ib/ft^,  a  choice  of  zero  or  "perfect”  bonding  between  the  pay- 
load  and  the  crushable  material,  and  a  choice  of  a  selected  balsa- like  or 
honeycomb- like  class  of  crushable  material.  Overall  radii  resulting  from  the 
designs  range  from  1.81  ft  to  about  3.45  ft,  and  the  resulting  crushable 
material  weights  vary  from  a  little  over  124  lb  to  almost  763  lb.  The  fol¬ 
lowing  conclusions  are  drawn  from  the  examples: 

1.  The  simplified  model  has  the  lower  crushable  material  weight  of  the 
two  models,  the  greatest  difference  being  approximately  15  percent 
for  an  example  without  payload  penetration. 
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2.  The  greatest  difference  between  the  two  models  is  a  29  percent 
discrepancy  in  the  maximum  crushing  stress  determined  by  a  design 
example  without  penetration  for  the  honeycomb-* like  class  of  material. 

3.  In  the  absence  of  payload  penetration,  the  crushable  material  weight 
and  the  overall  radius  increase  with  increasing  payload  radius,  that 
is,  with  decreasing  payload  packaging  density. 

4.  For  the  honeycomb -like  class  of  material  and  the  100- lb  payload  (with 
the  maximum  g  loading  held  at  2000  to  minimize  weight)  payload 
penetration  does  not  exist  for  the  detailed  model.  For  the  simpli¬ 
fied  model,  however,  penetration  slightly  reduces  both  the  crushable 
material  weight  and  the  overall  radius  at  the  lowest  payload  radius, 
0.6  ft,  the  only  radius  among  those  tried  for  which  penetration 
exists  with  honeycomb.  This  advantage  of  penetration  means  that 
bonding  should  be  avoided  with  the  present  honeycomb- like  material, 
for  which  the  specific  energy  absorption  increases  with  crushing 
stress  (a  variation  more  common  even  for  balsa  than  that  for  the 
present  balsa-like  class). 

5.  For  the  balsa- like  class  of  material  (with  a  payload  weight  of 
450  lb),  penetration  drastically  increases  the  crushable  material 
weight  at  the  lowest  payload  radius  (1.0  ft)  and  slightly  increases 

it  at  the  next  lowest  radius  (1.2  ft).  This  combines  with  nonpenetra¬ 
tion  results  to  give  a  minimum-weight  radius  for  an  unbonded  payload 
between  1.2  and  1.6  ft  and  a  corresponding  minimum-weight  payload 
packaging  density  between  26.33  and  62,15  Ib/ft^. 

6.  An  impressive  benefit  of  payload  penetration  is  a  decrease  in 
crushable  material  weight  by  a  factor  greater  than  4,  which  occurs 
for  the  100-lb  payload  with  a  radius  of  0.7  ft  when  the  (selected) 
honeycomb -like  class  of  material  is  required  without  penetration  (by 
the  g  loading  ceiling  of  2000),  but  when  the  more  efficient  balsa¬ 
like  class  is  feasible  with  penetration. 
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APPENDIX  A 


DEVELOPMENT  OF  BASIC  EQUATIONS 
GENERAL  GOVERNING  EQUATIONS 


Figure  2  shows  a  fairly  general  landing  geometry  at  the  start  of  impact 
and  at  two  later  instants  of  time  during  impact.  The  only  shape  requirement 
is  symmetry  about  two  mutually  perpendicular  vertical  planes,  and  this 
requirement  is  made  for  compatibility  with  the  assumption  of  pure  vertical 
translation.  The  crushing  strength  and  density  of  the  crushable  material  are 
also  required  to  be  symmetrical  but  otherwise  generally  variable.  The  land¬ 
ing  surface  is  assumed  to  be  perfectly  flat  and  rigid.  The  payload  is 
assumed  to  be  perfectly  rigid  and  is  considered  unbonded  to  the  crushable 
material  in  determining  the  start  of  relative  motion  (if  any) ,  herein 
referred  to  as  "payload  penetration." 

The  analysis  begins  with  the  phase  of  impact  shown  in  figure  2(b).  In 
this  phase,  the  crushable  material  has  begun  to  crush  against  the  landing 
surface,  but  payload  penetration  has  not  begun.  Hence, 


qp  =  q 


(Al) 


where  qp  and  q  are  the  absolute  displacements  of  the  payload  and  the 
un crushed  crushable  material,  respectively,  as  shown  in  figure  2(b).  The 
equation  of  motion  for  the  constant  payload  mass  mp^  is  written  as: 


mq  =  mg-  /  jo  dA  r\o'\ 

po^  po^  J  J  Aqpo  (A2) 

^po 

where  q  is  the  second  time  derivative  of  q,  g  the  local  acceleration  due 
to  gravity,  cjypo  the  vertical  component  of  normal  stress  on  a^nd  Apo 

the  area  over  which  the  stress  acts.^ 


The  equation  of  motion  is  now  written  for  the  variable  mass  m^o  ’ 
of  uncrushed  crushable  material  in  figure  2(b)  under  the  assumption  that  each 
particle  in  m^o  "  ^cl  is  moving  at  the  same  vertical  velocity  q  at  a 
given  instant.  The  equation  is 

^Note  that  the  vertical  component  of  shear  stress  on  the  payload  could 
have  been  added  to  cxypo  equation  (A2)  and  in  the  later  equations  con¬ 

taining  Uypo.  There  are,  however,  no  Uypo  terms  in  the  major  governing 
equations  except  for  the  equation  that  determines  the  existence  of  payload 
penetration.  When  penetration  is  about  to  start,  the  most  important  stresses 
are  cTypo  values  that  are  almost  large  enough  to  crush  the  material  (in  its 
strongest  direction  if  it  has  one  and  has  been  deployed  to  utilize  the  fact) . 
Shear  stresses  are  negligible  by  comparison  and  are  often  incorporated  in  the 
tests  to  determine  the  crushing  strength. 
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where  Oy  is  the  vertical  component  of  the  "mostly  static"  and  "mostly 
normal"  stress  capable  of  deforming  the  crushable  material  plastically  oyer 
the  area  term  "mostly  static"  implies  the  assumptions  that  the 

damping  force  is  incorporated  and  that  the  coupling  between  vertical  and 
horizontal  velocity,  or  dynamic  buckling  effect^  is  also  incorporated.  The 
term  "mostly  normal"  implies  the  incorporation  of  a  small  shear  stress.  (See 
the  discussion  of  figs.  1(b)  and  1(c)  imder  "Properties  of  Typical  Crushable 
Material.")  The  area  Aci  is  the  intersection  between  -  mci  and  the 

mass  of  crushed  material  m^^  in  figure  2(b).  The  crushed  material  is 
assumed  to  have  been  compacted  sufficiently  to  deliver  the  required  crushing 
stress  Oy;  and  successive  layers  of  material  are  transferred  from  m^o  " 
to  m^i  at  the  boundary  A^i  as  they  undergo  plastic  deformation.  (Such  a 
process  is  observed  experimentally  in  the  crushing  of  plastic  foam,  balsa 
parallel  to  its  grain,  and  honeycomb  parallel  to  its  axis.) 


The  assumption  of  a  uniform  (though  time  variable)  vertical  velocity  q 
throughout  m^o  •"  two  implications  besides  the  crushing  by  layers 

just  described.  First,  the  elastic  stress  waves  that  establish  a  uniform  q 
must  obviously  travel  at  speeds  far  greater  than  q  (i.e,,  q  must  be  a  low 
subsonic  value) .  Second,  the  deformation  waves  associated  with  the  elastic 
stress  waves  must  be  small  enough  not  to  affect  the  magnitude  of  q. 

Equations  of  motion  could  also  be  written  for  m^i  in  figure  2  and  for 
the  elemental  layer  of  mass  being  transferred  to  mci*  These  equations 
would  determine  dynamic  stresses  in  the  sense  of  variable  mass,  however,  and 
are  not  needed  since  Oy  is  assumed  to  be  a  measured  quantity. 

Equations  (A2)  and  (A3)  can  be  used  to  determine  whether  payload 
penetration  occurs  and,  if  so,  at  what  displacement.  If  both  equations  are 
solved  for  q  -  g  and  the  results  equated  and  rearranged,  the  result  is: 


where  +  m^Q  -  is  the  total  time-dependent  mass  at  velocity  q  in 

the  absence  of  payload  penetration.  Equation  (A4)  is  the  only  major  governing 
equation  containing  Uypo;  the  neglect  or  incorporation  of  shear  stresses  is 
justified  in  footnote  1  of  this  appendix. 

If  equation  (A4)  is  considered  formally  solved  for  crypo^  payload 
penetration  cannot  begin  until  those  stresses  become  large  enough  to  cause 
plastic  failure,  that  is,  become  a  Oy  distribution.  This  will  occur  only 
if  the  geometry,  masses,  and  impact  conditions  are  such  that  the  displacement 
q  becomes  large  enough  to  bring  A^-^  up  to  the  necessary  size.  It  can  be 
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seen  in  equation  (A4)  that  the  following  quantities  tend  to  prevent  or  post¬ 
pone  payload  penetration:  a  large  payload  bearing  area  Apo,  a  small  payload 
mass  mpo,  a  large  total  mass  mpo  +  m^o  -  iUci,  and  a  small  bearing  area  of 
compacted  material  A^.^  • 

Equations  (A2)  and  (A3)  can  be  added  to  give  the  major  governing 
equation: 


(m  +  m  -mid 
po  CO  cl  ^ 


(m  +  m  -  m  , )  g 

po  CO  cl  * 


// 


a  dA 

V 


(A5) 


‘Cl 


This  is  the  simplest  equation  of  motion  to  apply  during  the  impact  period 
prior  to  payload  penetration.  It  does  not  contain  Oypo  (see  footnote  1  of 
this  appendix) . 

If  equation  (A4)  shows  that  payload  penetration  has  begun,  the  geometry 
of  figure  2(c)  applies.  The  constant  payload  mass  nipo  is  then  positioned 
by  the  coordinate  qp  so  that  equation  (A2)  becomes: 


m  q 
po  P 


m  g 
po* 


-// 


a  dA 
vpo 


(A6) 


^po 


The  mass  mpi  in  figure  2(c)  consists  of  compacted  crushable  material. 
The  velocity  qp  of  mpo  is  assumed  to  apply  uniformly  throughout  mpi . 

Hence,  the  equation  of  motion  for  mp]^  can  be  written: 


m 


pi 


m  g  + 

pi 


//°vpo"-//-v' 


dA 


(A7) 


po 


pi 


where  Ap^  is  the  intersection  between  mpi  and  the  variable  mass 
^co  ”  ^cl  "  ^pi  uncrushed  crushable  material,  and  Oyl  differs  from  the 
"mostly  static"  crushing  vertical  component  Oy  acting  over  Ap^  only 
because  the  elemental  layer  of  mass  dmpj  is  transferred  to  ii^pi  • 


given 


The  equation  of  motion  for  dmp^  becomes,  with  the  definitions  just 


dm 


qp  q 


plV  dt 


-)■/  /%'“-/  /%“ 


or 


Cqp  -  q)mpi  =  /  I  -  J  I 


(A8) 
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where  nipj  =  dmpi/dt,  where  the  higher-order  quantity  gdirip^  has  been 
neglected,  and  where  qp  and  q  represent  the  final  and  initial  velocities 
of  dmpi . 

Equations  (A6) ,  (A7) ,  and  (A8)  can  be  added  to  give 


Api 


(A9) 


A  second  major  equation  of  motion  defines  the  forces  on  the  variable  mass 
i^co  “  "  ^pi  of  uncrushed  crushable  material  in  figure  2(c).  Under  the 

assumption  that  each  particle  in  m^o  -  ^cl  "  moving  at  the  same 

vertical  velocity  q,  this  equation  is 


-  "'cl  -  =  Ko  -  “cl  - 


where  Api  is  defined  as  in  equations  (A7) 
analogously  to  that  in  equation  CA3) . 

Equations  (A9)  and  (AlO)  define  the  motion  during  payload  penetration 
since  all  quantities  can  be  defined  in  terms  of  q  and  qp  and  their  deriva¬ 
tives.  The  two  major  governing  equations  are  coupled  in  the  most  general 
case  by  the  qp  -  q  term,  the  dependence  of  mpi  on  qp  -  q,  and  the 
dependence  of  and  of  Oy  over  Ap^  on  qp  -  q.  Equations  (A9)  and 

(AlO)  do  not  contain 


ih-fj 


a  dA 

V 


Aci 


(AlO) 


and  (A8)  and  A^i  is  defined 


GOVERNING  EQUATIONS  FOR  ZERO  SHEAR  DEFORMATION, 
UNIFORM  COMPACTING  STRAIN,  AND  UNIFORM 
MATERIAL  DENSITY 


Although  figure  2  does  not  show  any  shear  deformation,  the  equations 
derived  so  far  would  apply  even  if  such  deformation  were  present.  The  speci¬ 
fication  of  Oy,  however,  requires  knowledge  of  the  location  of  the  surfaces 
Aci  and  Api,  and  shear  deformation  is  ruled  out  at  this  point  to  retain  the 
simple  surfaces  implied  by  figure  2.  Then  Aci  ^pi  determined  by  the 

height  variables  hci  and  hp^ ,  where  h^^  is  the  local  height  of  m^i  in 

figures  2(b)  and  2(c)  and  where  hp^  is  the  local  height  of  mpi  in  fig¬ 
ure  2(c).  If  shear  deformation  had  been  considered,  there  would  be  a  trail¬ 
ing  of  material  to  make  mpi  wider  than  mpo  and  a  possible  lifting  off 
the  ground  of  the  edges  of  m^i. 
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The  actual  lengths  of  h^i  and  h^.^  are  evaluated  in  terms  of  a  readily 
measurable  variable  e.  This  variable  is  defined  as  the  compressive  compact¬ 
ing  strain  of  the  crushable  material  ('^compacting  strain"  is  the  strain  at 
which  the  crushing  stress  rises  abruptly  from  a  relatively  constant  value  for 
a  test  specimen  of  uniform  cross  section) .  The  effect  of  the  direction  of 
maximum  strength  is  ignored  under  the  assumption  that  the  material  will  com¬ 
press  to  the  same  compacted  strain  along  any  axis. 


The  use  of  e  is  made  feasible  by  the  assumption  that  the  energy 
absorbing  process  consists  of  the  vertical  crushing  of  separate  vertical  rods 
of  material.  Then  h^  in  figures  2(b)  and  2(c)  is  the  total  shortening 
deformation  of  a  rod  having  an  initial  length  of  he  +  hci,  and  qp  q  is 
the  corresponding  deformation  of  a  rod  having  an  initial  length  of 
qp  -  q  +  hpi,  where  hpi  is  the  local  height  of  in  figure  2(c). 

Since  the  variable  e  di  is  (by  the  definition  of  strain)  the  total  shorten¬ 
ing  deformation  of  each  successive  rod  element  dil  to  be  crushed,  it  is  seen 
that 


qp  -  q 


^c 


qp-q+hpi 

/ 


/ 


h  +  h  , 
c  cl 


z  dl 


(All) 


Under  the  assumption  of  uniform  compacting  strain  (e  =  =  constant), 

the  equations  just  derived  can  easily  be  solved  for  hp^  and  h^-j 


1  -  £1 


A 


^pl 


cl 


TiT  '’p  - 


1  “  e. 


(A12) 


The  effects  of  equations  (A12)  for  the  constant  can  be  seen  in  figure  3, 

which  is  otherwise  identical  to  figure  2.  In  figures  3(b)  and  3(c),  the  mass 
m^i  has  become  a  foreshortened  image  of  the  volume  V^,  which  would  lie 
beneath  the  landing  surface  if  there  had  been  no  crushing;  and  the  height 
hpi  of  the  mass  has  become  constant  in  figure  3(c). 


The  next  limitation  to  be  imposed  in  the  analysis  is  that  of  uniform 
density  of  the  crushable  material  prior  to  crushing.  Then  m-pi  and  m^  in 
equations  (A4) ,  (A5) ,  (A9) ,  and  (AlO)  become: 
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m 


cl 


Pck 


+  V  ) 
P 


Pck 


+  V  ) 
c 


where  Vpj  is  the  volume  of  nipi ,  Vp  the  volume  swept  out  by  mpo  during 
payload  penetration,  and  Vci  the  volume  of  m^i .  When  Vp^^  is  related  to 
Vp  and  Vci  to  according  to  equations  (A12) ,  the  masses  mpj  and  mci 
become 


1  -  £1 


m 


ck 


+  V 


j 


(A13) 


The  relations  between  the  volumes  Vp  and  V^  and  the  displacements 
and  qp  can  be  determined  by  the  following  volume  formulas: 


V  =  A  ,  (q  -  q) 
p  poh^^p 

0  J 


q 


(A14) 


where  Ap^j^  is  the  horizontal  planar  projection  of  Apo  and  Apj ,  and  Ash  is 
a  horizontal  cross  section  in  V^. 


Another  useful  set  of  relations  involves  the  stresses.  The  stress 
integrals  in  the  major  governing  equations  (A4) ,  (AS),  (A9) ,  and  (AlO)  can  be 
written 

/  /°vpo  "  -f  j%o  " 

A  A  . 

po  poh 


J 


(A15) 
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where  Opo  is  the  normal  stress  on  the  payload  prior  to  penetration,  a  is 
the  normal  crushing  stress  on  any  surface,  and  Acih  is  the  horizontal 
planar  projection  of  Aci .  For  integration  of  the  right-hand  sides  of  equa¬ 
tions  (A15) ,  the  normal  stresses  must  be  evaluated  at  the  intersections 
between  the  curved  surfaces  (Apo,  Ad,  Api)  and  vertical  lines  through  the 
centroids  of  the  elements  dA  in  the  horizontal  planar  surfaces  (Apoh,  Adh, 
^poh)  • 

A  final  useful  relation  is  the  following  change  of  variable: 


With  equation  (A16) 


“  ^  ^  ^  dU  _  dU  _  1  dCU^j 
^  dt  “  dt  dq  2  dq 


(A16) 


(A17) 


Equations  (A13j  for  uniform  density  and  uniform  compacting  strain, 
together  with  the  volume  equations  (A14)  and  the  stress  equations  (A15) ,  are 
now  introduced  into  the  basic  governing  equations.  In  addition,  one  of  the 
governing  equations  is  reduced  to  first  order  by  equation  CA17).  Thus  equa¬ 
tion  (A4),  which  determines  the  start  of  payload  penetration  (if  any), 
becomes : 


If 


m 


po 


If 


a  dA 


0  dA  = 
po 


clh 


poh 


m  +  m  -  (p  1  /e,  ) 
po  CO  ^^ck'  k-^ 
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A  ,  dh 
sh  c 


Equation  (A5) ,  which  defines  the  impact  prior  to  penetration,  becomes 

I 


1  d(U2) 

2  g  - 


clh 


m  +  m  -  (p  1  /gi  ) 
po  CO  ^  ck  k 
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A  ,  dh 
sh  c 


(A18) 


(A19) 


Finally,  equations  (A9)  and  (AlO)  for  payload  penetration  become; 
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_ p°h 

“po  "  ‘'■ck^poh^'k’  %  ■ 


q  =  g  - 


/  fa  dk  -  f  I  a  dk 

\oh 


(A20) 


cih 


m 


CO 


Equations  (A18]  through  (.A20)  define  the  problem  for  the  landing  geometry  of 
figure  3,  providing  the  stress  and  volume  integrals  can  be  evaluated. 


EVALUATION  OF  VOLUME  AND  STRESS  INTEGRALS 
FOR  SPHERICAL  GEOMETRY 


The  evaluation  of  the  volume  and  stress  integrals  is  facilitated  by 
restriction  of  the  landing  configuration  to  a  simple  shape.  The  sphere  is 
particularly  useful  for  landing  packages  and  crushable  coverings  that  must 
absorb  energy  from  impacts  in  any  direction.  Hence,  the  geometry  is  now 
specialized  to  concentric  spheres. 

The  spherical  landing  package  is  shown  in  figure  4(a)  at  the  start  of 
impact,  in  figure  4(b)  during  impact  but  prior  to  payload  penetration,  and  in 
figure  4(c)  after  penetration.  All  angles  and  radial  lengths  are  shown  in  a 
plane  of  maximum  value. 

The  volume  integral  in  equations  (A18)  through  (A20)  can  be  evaluated 
immediately  on  the  basis  of  sketch  (a),  which  is  a  simplification  of  figure 
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4(a)  or  4(b).  Thus 


or 


Vq  =  ^  dh^  =  J*  Trr^  ^  f  " 

0  0  *^0 

''c^/''*sh  <ihc  -  q) 

0 


(R  -  q  +  h^)2]cih^ 


CA21) 


It  is  also  noted,  for  certain  volume  terms  and  integration  areas  in  equa¬ 
tions  (A18)  and  (A20) ,  that 

Apoh  ==  TiRp^  (A22) 

as  can  be  seen  from  figure  4  and  the  definition  of  Apoh* 

One  of  the  stress  integrals  can  be  evaluated  immediately  under  the 
temporary  assiomption  that  Opo  =  apok  =  constant.  Then,  with  equation 
(A22)  , 

/I  0  dA  =  a  1 A  ,  =  ttR  ^0  ,  (A23) 

J  po  pok  poh  p  pok 

^poh 


The  assumption  that  Cpo  =  cfpok  i^^y  be  valid  only  when  Opok  approaches 
the  normal  crushing  stress  a,  that  is,  when  payload  penetration  is 
approached.  Fortunately,  this  is  the  only  region  of  interest  for 
equations  (A18)  and  (A23) . 

Prior  to  evaluation  of  the  stress  integrals  for  a,  another  limitation 
is  applied  to  the  analysis,  namely, 

0  =  00  COS  a  ;  a  <  90°  (A24) 

where  a  is  the  angle  between  the  local  normal  to  the  stressed  area  and  a 
radial  line  from  the  center  of  the  spherical  system  (as  if  un deformed)  through 
the  point  of  stress  application,  (The  rationale  for  this  limitation  is  given 
immediately  prior  to  eq.  (2).) 

The  next  step  is  to  determine  a  for  mci  in  figure  4(b)  (with 
applicability  after  payload  penetration  as  well).  Any  point  on  the  surface 
of  revolution  shown  by  the  dotted  portion  of  the  circle  is  defined,  with  the 
aid  of  sketch  (a) ,  by 


r^  +  (R  -  q  +  hc)^  = 
From  the  second  of  equations  (A12) 
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and  h^i  can  be  determined  from  equation  (A25) .  Then,  after  some 
manipulation,  together  with 


// 


clh 


we  obtain 
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a  dA  =  I  a27Tr  dr 


2'rT[(R  -  q) -  (1  -  e.)R^]r  dr 


(1  - 


(r2  -  r2)]|r2  +  (1/e^2)[r  .  q  -  (1  -  e^)/R2-  ^ 


/2 


With  the  substitution  b  -  r^,  this  becomes 


If 


[(R  -  q)b  -  (1  -  eR)R2]b  db 


a  dA  =  27Ta 


clh 


:-,(l 


-  (1  - 


)2R2jjR2  -  b2  +  (1/£^2)[R  -  q  -  (1  -  C^)b]2|^ 


1/2 


(A27) 

and  equation  (A27)  is  the  stress  integral  for  the  upper  surface  A^i  of 
mci  in  figure  4(b)  or  4(c).  The  integral  contains  no  singularities  and 
hence  can  readily  be  evaluated  numerically. 


The  angle  a  is  now  determined  for  m-pi  in  figure  4(c).  Thus,  with 

sketch  (c)  for  detail. 
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or,  with  the  first  of  equations  (A12) , 


a  =  arc  sin  -r - arc  tan 

K-n 


(i/ev)(qr.  -  q)  + 


(A28) 


When  equation  (A28)  is  introduced  into  equation  (A24) ,  the  normal 
crushing  stress  a  is  determined.  When  a  is  put  into  o  dA  with 


R, 


Apoh 


the  result  is 


/p 


If 


a  dA  =  a 


poh 


2»[(qp  -  -  rp2  *  Rp^jrp  top 

Rp/2[(qp  -  q)/e|^]/Rp2  -  Tp^  .  ( (qp  -  q)/e^]2  t  R 


With  the  substitution  bp  =  /kp^  -  rp^,  this  becomes 


r  r  2TTa  i 

A  -  J 

poh  0 


{[(qp  -  q)/Ek]^  *  % 

/2[(q  -  q)/ej^]b  +  [  (q  -  q)/£j^]2  +  R 


(A29) 


Equation  (A29)  is  the  stress  integral  for  the  lower  surface  Api  of  mpi  in 
figure  4(c).  It  can  be  integrated  directly  to  give 


// 

%oh 


2iTa 


a  dA  = 


15R, 


Rp^  *  3Rp“ 


(^) 


-v](. . 


(A30) 


Since  equation  (A30)  is  indeterminant  at  [(qp  -  q)/e]  =  0,  that  is,  at 
the  start  of  payload  penetration,  the  best  way  to  evaluate  ^^^*0  dA  is  to 

integrate  equation  (A29)  numerically.  At  [(qp  -  q)/e]  =  0>  equation  (A29) 
becomes 
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If 


=  ttR 


(A31) 


[(qp-q)/^]=o 


Comparison  of  equation  (A31)  with  equations  (A23)  and  (A24)  serves  as  a 
partial  check  on  the  development  of  equations  (A29)  through  (A31)  since 
a  =  0°  (or  a  =  Oq)  on  the  payload  at  the  start  of  penetration. 


Dimensionless  Governing  Equations  for  Spherical  Geometry, 

With  Termination  Conditions 

Governing  equations  that  have  been  derived  in  this  appendix  are  now 
made  dimensionless  in  a  form  convenient  for  numerical  solution  of  the  impact 
problem  with  spherical  geometry.  Prior  to  payload  penetration,  equation 
(A19)  applies  and  becomes,  with  equations  (A21)  and  (A27) , 


_  1  d(w^)  __  g 


Fcl(z) 


(A32) 


where  ge  is  the  acceleration  due  to  gravity  on  earth,  npge  and  nn,(jge  are 
the  actual  and  maximum  design  payload  decelerations,  respectively,  and 


i^mdSe^ 


z  =  a. 

^  "  R 


(A33) 


Kn  =  - 
K  m 


po^idS( 


The  dimensionless  total  mass  m(z)  in  figure  4(b)  is  defined  by 


m(z) 


^po  ^CO  “  ^cl 


=  1  + 


^po  \^pR 


z^(3  -  z) 


(A34) 


(4/3)  ttR^ 


(A35) 
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(A37) 

Equations  (A32)  and  (A36)  apply  until  z  reaches  a  value  Z3,  at  which 
payload  penetration  (if  any)  starts.  The  start  of  penetration  is  defined  by 
equation  (A18)  with  equations  (A23) ,  (A34) ,  and  (A36) .  Thus, 


R  Y 

Rp/  m(z) 


CA38) 


and  penetration  starts  (z  =  Zs)  if  and  when  Opok  reaches  Oq  (providing  the 
payload  is  not  attached  to  the  crushable  material) . 

At  the  start  of  penetration,  it  becomes  convenient  to  introduce  three 
new  variables  and  one  constant  as  follows: 
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_  q  R 


X  E  t 


md^e 


qp  -  q 


(A39) 


^Rp^Oo 


^  '"poVd^e 

The  initial  conditions  for  penetration  are 


y(x  =  0)  =  Yg  E  Zg  ^ 


dx  )  I  dx 


e(x  =0)  =  e„  =  0 


(A40) 


I  M  ° 

x=0  \  /s 


where  Wg  is  the  value  of  w  determined  by  equation  (A37)  at  z  =  Zg. 

During  penetration,  equations  (A20)  apply,  with  equations  CA21) ,  (A22) , 
(A29),  (A36),  and  (A39) .  Thus 


d^e  ^  ^  ^  g 

j  2  ,2  n  jg 

dx  dx  md^e 


(TTRp3pcke]^/mpo)  (de/dx)^  +  KpEp^Ce) 
mpenCe) 


d^y  _  g  ■  ^p^po^®^ 

dx^  ^mdSe  m„„(y,e) 


(A41) 
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where  Fci (y)  is  determined  by  equation  (A36)  with  (Rp/R)y  substituted  for 
z;  mpej^(e)  is  the  dimensionless  total  penetrating  mass  in  figure  4Cc)  defined 

by 


^po  ^pl 


3r  ck 


1  +  ^  —  e 

“Vp  / 


(A42) 


"  3 

(4/3)TrRp" 

(y,e)  is  the  dimensionless  crushable  mass  in  figure  4(c)  defined  by 


(A43) 


m  (y.e)  = 

pr'-'  ’ 


‘CO  -  “ci  -  >1  "*co  1  2A  L.  - 

"•po  ’  4ek  \Pp  /  V  / 


(A44) 


and  Fpo(e)  is  the  dimensionless  crushing  force  on  mpen(yp,y)  defined  by 


%R  J 

P  O 


(esp  +  l)Sp  dSp 


2esp  +  e  +  1 


(A45) 


with  Sp  =  bp/Rp. 

Equations  (A32) ,  (A37) ,  (A40) ,  and  (A41)  are  the  required  governing 
equations  and  initial  conditions.  With  rebound  excluded  and  with  no  payload 
penetration,  the  impact  is  terminated  (z  =  z^ax)  when  w  =  0.  Thus 


w(zmax)  =  0 


(A46) 


With  rebound  excluded  but  payload  penetration  present,  the  impact  is  termi¬ 
nated  (y  =  ymax  ^nd  y  +  e^e  =  (y  +  £k®)max^  when  dy/dx  =  0  and 
dy/dx  +  ek(de/dx)  =  0.  Since  e  increases  monotonically  with  time  for  the 
present  examples,  the  termination  condition  becomes 


y=ymax 


(A47) 


As  the  payload  begins  to  penetrate  the  crushable  material,  equations 
(A41)  are  integrated  simultaneously.  This  continues  until  the  first  of 
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equations  (A47)  is  satisfied,  thereby  providing  initial  conditions  for  the  next 
phase  of  the  problem.  In  this  next  phase,  the  second  of  equations  (A41)  can 
be  ignored  and  the  first  equation  solved  alone.  The  uncoupled  solution  of 
the  first  of  equations  (A41)  continues  until  the  second  of  equations  CA47)  is 
satisfied. 


Simplifications  for  Constant  Mass,  Infinitely  Thin  Crushed  Material, 
and/or  Constant  Penetration  Resistance 

The  governing  equations  are  now  written  with  constant  mass  (p^k  =  0) 
Equations  (A32)  and  (A34)  are  combined  to  give 


^  d{w^)  _  g 
^  "  2  dz  H”? 


FciCz) 


(A48) 


where  Fd  (z)  is  given  by  equation  (A36) ,  Equation  (A37)  becomes 


z  z 

W2  =  Wo2  +  2J  dj  =W^2+2y* 


%dge  "  1  +  (mco/™no) 


=  Wo  + 


2gz 

^md®e 


2Kr  r 

1  +  (m  /m  )  / 
CO  po^  J 


(A49) 


and  equation  (A38)  gives 


^Vok  / R  V  Pci  (2;) 

\  P/  po 


(A50) 


After  payload  penetration  starts,  equations  (A41)  apply,  and  p^k  =  0  gives 


d^e  d^ 


md  dx  dx  md^e  ^ 

dV  _  g  -  KpFpo(e) 

"md^e  "'co/™po 


(A51) 


where  F(>i  (y)  is  given  by  equation  (A36)  with  y(Rp/R)  substituted  for  z  and 
where  FpQ(e)  is  given  by  equation  (A45) . 

Regardless  of  whether  pck  =  0  or  0,  equation  (A36)  is  greatly 

simplified  by  the  specialization  to  er  =  1  Cidfinitely  thin  crushed  material) 
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Thus,  for  E]^  =  1,  equation  (A36)  can  be  integrated;  and  z  =  y(Rp/R)  can  be 
substituted  as  follows: 


Fci(z)  =  2z(l  -  z) 


(A52) 


Regardless  of  whether  pck  =  0  or  0  and  whether  =  1  or  5^  1,  equation 
(A45)  is  greatly  simplified  by  the  assumption  that  the  payload  stress  and 
force  maintain  their  initial  penetration  values  throughout  the  penetration 
stroke,  that  is,  the  assumption  that  Fpo(e)  can  be  represented  by  FpQ(O), 
which  is  unity.  Thus,  the  assumption  has  the  form 

Fpo(e)  =  Fpo(O)  =  1  (A53) 

Equations  (A52)  and  CA53)  can  be  used  individually  or  together  in  the  govern¬ 
ing  equations  for  p^k  0  or  for  p^k  =  ^  whenever  the  simplifying 
assumptions  are  appropriate. 

If  the  assumption  of  equation  (A53)  is  made  (constant  penetration 
resistance)  and  if  Pck  0  (constant  mass),  equations  (A41)  can  be  decoupled 
by  a  change  of  variable,  regardless  of  whether  ek  =  1  or  ^  1,  as  can  equa¬ 
tions  (A51) ,  which  are  specialized  for  Pck  =  0-  '^^e  governing  equations 

(A48)  through  (A51)  will  be  rewritten,  however,  for  ek  =  1,  that  is,  equa¬ 
tions  (A52) ,  as  well  as  the  assumption  of  equation  (A53) .  The  change  of 
variables  is 
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(A55) 


Cy  -  Ys) 
(A61) 

where  Vs  is  the  dimensionless  velocity  corresponding  to  y^. 

Equations  (ASS)  through  (A61)  are  the  governing  equations  for  infinitely 
thin  crushed  material  (ej^  =  1)  and  constant  mass  (p^.)^  =  0),  with  the  assump¬ 
tion  of  constant  penetration  resistance  (Fpo(e)  =  1)  required  after  payload 


+  2  - 

b 


m  /m 
CO  po 


R  ,  .2 

f>'  *  >'s>  *  3 
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penetration  but  not  before.  In  the  absence  of  penetration,  equation  (A56) 
can  be  used  to  determine  y^ax  setting  v  =  0.  Thus  Ymax  solu¬ 

tion  of  the  cubic  equation 


1  +  (mco/mpo) 


^md^e 


0 


(A62) 


The  dimensionless  velocity  Vg  at  which  penetration  starts  (if  any)  can  be 
determined  by  writing  equation  (A56)  as 


2gys 


2Kpys2[(R/Rp)  -  (2/3)y3] 

1  +  (m  Tin  ) 

^  CO  po-^ 


(A63) 


where  is  the  solution  of  equation  (A58) .  Then  the  maximum  payload 

dimensionless  displacement  yp  with  penetration  is  determined  by  equa¬ 
tions  (A60)  and  (A63)  with  Yp'^^^O  as 


V  2  .. 


yp 


max 


2Kpys(l  -  jys[CR/V  ■  (2/3)ys]/[l  +("100/^31) 

2[Kp  -  (g/nn,dge)J 


(A64) 


A  similar  procedure  with  equations  (A61)  and  (A63)  and  with  v  =  0  gives 
ymax  as  the  solution  of  the  cubic  equation. 


2K- 


^  ^  (y^  +  y  y  +  y  ^ 


m  /m  13  ^-^max  -^max^s  ^s 
CO  po  [ 


^  ^  ^^^max  ^  ^s^  ^^max  " 


(y  “  y  )  ■’■ 


2gy 


max 


+  V, 


2Kpy/[(R/Rp)  -  (2/3)y^] 


max  S'  n  ,g 
md^’e 


1  +  (m  /m  ) 

CO  po 


=  0 


(A65) 

It  should  be  noted  that  yp^ax  generally  sufficient  for  design  problems 
involving  payload  penetration  and  that  equation  (A65)  for  ymax  included 
only  for  completeness  and  possible  checking  purposes. 
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APPENDIX  B 


SAMPLE  CALCULATION  FOR  DESIGN  BY  SIMPLIFIED  MODEL 
WITHOUT  PAYLOAD  PENETRATION 


Under  the  assumptions  that  Pck  =  0,  =  I,  and  g  =  0  (simplified 

model) ,  a  calculation  is  performed  in  this  appendix  to  yield  a  design  of  a 
crushable  casing  for  the  following  conditions: 


Uq  =  300  ft/sec 
ed  =  0.7 

Em  =  0.8 

%ax  ■ 

Wpo  =  nipoSe  ^^O  lb 
Rp  “0.6  ft 

Payload  penetration  prevented  by  bonding 
if  necessary 

The  calculation  corresponds  to  case  1  in  table  1,  as  indicated  by  the  last 
four  of  the  above  conditions.  The  calculation  follows: 


°  _  (300)2 

Cnpmaxge)Rp  (2000) (32.17) (0.6) 


2.331 


From  figure  8  for  GjJ  =  0.7 


0.286 


or 


0.6 

0.286 


2.097  ft 


From  figure  9  for  ed  =  0.7  with  Rp/R  =  0.286, 

■^mo  “  2.74 
Uo" 

check:  y- - —75-=  (2.331)  (0.286)  =  0.6665 

^"Pmax^eJR 


From  figure  6 


^max 


"max 

R 


0.500 
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or 


q  =  0.500(2.097)  =  1.049  ft 
^max 

From  figure  7  for  =  0.7 

Jma  =  2.74  (check  is  perfect  to  3  places;  therefore  use  2.74) 
From  figure  12  with  Rp/R  =  0.286, 

^ma  ”  55.5 

^mJma  =  0.8(2,74)  =  2.192 

=  2.503«10--  ft-> 

^0  (300) 


•^ma^'^po 


N^,.Rp3(144) 

From  equation  (27)  with  in  psi 


(2.74)  (100)  _  r.c;Wft 

(55.5)it(0.216)  (144)  0.5053  10  psi/ft 


2  1 Q2 

SEA  =  - -  ft-lb/lb 

2.503x10-'^  -  (0.5053x10-1/0 


Now  decide  which  material  to  use  in  figure  5  (trial  eliminates  the  balsa- like 
class,  which  is  too  strong  for  low  payload  weight).  From  figure  5  for 
honeycomb- like  material, 

try  SEA  =  10,400  ft-lb/lb  at  Oq  =  1,000  psi 

2  192x10^ 

°  2.503  0.5053  ° 

try  SEA  =  10,870  ft-lb/lb  at  =  1,100  psi 
SEA  =  _ia92><10^^  jQ^720  ft-lb/lb 


2. 

503 

-  0. 

5053 

10 

,870 

ft- 

lb /lb 

2.192x10 

2. 

503 

-  0. 

459 

10 

,800 

ft- 

lb /lb 

2.192x10 

4 

2.503  -  0.463 
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try  SEA  =  10,770  ft-lb/lb  at  =  1,080  psi 


•  2.5offo°l68  = 


therefore 


Oq  =  1,080  psi 
SEA  =  10,770  ft-lb/lb 

From  equation  (8),  with  Oq  in  psi  and  SEA  in  ft-lb/lb 

p  g  _  (0.8)  (1441  ri, 0801 

cm^e  SEA  10,770 


55  Ib/ft 


From  equation  (17) 


itipo  +  m^Q 


'JmaUo^^m 

ge(SEA)e 


(2.74)  (300)2(0.8) 
(32.17) (10,770) (0.7) 


0.8135 


m~o  0.8135 


-  1  =  1.229  -  1  =  0.229 


therefore 


4.366 


Wj,Q  =  4.366(100)  =  436.6  lb 
''^po  +  Wco  =  100  +  436.6  =  536.6  lb 


Check  by  recalculating  PcmSe-  From  equation  (25), 


N  3  (55.5)  (•^)  (0.216)  " 

(Check  is  adequate,  use  11,56  Ib/ft^.)  Check  to  see  if  bonding  is  necessary 
to  prevent  payload  penetration 


RpV  /  ^  \ 

Fy  y  (0.286)2(5.366)  =  0.4386 


From  figure  10 
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Zg  =  0.325 


^max  _  *^max  _  ‘^Pmax  _  0 .500 
Zs  %  ~  % 


1.538  >  1.000 


Therefore  bonding  is  necessary  to  prevent  payload  penetration. 
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APPENDIX  C 


SAMPLE  CALCULATION  FOR  DESIGN  BY  SIMPLIFIED  MODEL 
WITH  PAYLOAD  PENETRATION 


A  sample  calculation  is  given  for  an  approximate  design  of  a  crushable 
casing  under  the  assumptions  that  =  0,  =  1,  and  g  =  0  (simplified 

model)  together  with  FpQ(e)  =  1  for  penetration.  The  design  conditions  are 
the  following: 

Uq  =  300  ft/sec 
cd  =  0.7 
Em  =  0-8 

np  ^  2000 

i^max 

Op  =  1200  psi  for  balsa- like  material 

'^00  =  m  g  =  450  lb 
po®e 

Rp  =  1.2  ft 

Payload  penetration  permitted  (unbonded) 

The  last  four  conditions  indicate  that  case  11  of  table  1  is  being  calculated. 
From  figure  5  for  balsa-like  material. 


SEA  =  24,000  ft-lb/lb 


From  equation  (8),  with  Op  in  psi  and  SEA  in  ft-lb/lb 


p  g  =  =  10,8)  (144)  (1,200)  ^  3 

^cm^e  SEA  24,000 


From  equation  (22)  with  in  psi  (for  penetration) 


n^ 


_  7rRp2(144  Op) 


Pmax  m  g  450 

po^e 

For  iteration  with  figures  11  and  12,  calculate 


44) (144) (1,200) 


1,737 


1  /  \  1  (300)2 

Upn^ax^e^  r  (1737)  (32.17)  (1.2) 


0.959 
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1  +  0.959 


1.959 


(1  +  0.959)^  =  3.838 


SEA 


24,000 


nr,  e  R  1737(0.8)  (1.2) 

Rmax  m  p  j  k  j 

">00 

first  - =  1.0  (guess) 

m 

po 


=  14.385 


first  N 


=  1 
mu  3.838 


0.5213 


From  figure  11  for  =  0.7 


From  figure  12 


first  =  0.8435 

..  ^  0.8435  _ 

R-= 


first  Njjig  =  15.5 


m 


second 


CO  15,5 


m  14.385 
po 


1.078 


second  N 


mu 


1  +  1.078 
3.838 


=  0.5415 


From  figure  11  for  =  0,7 


From  figure  12 


second  =  0.838 

,  0.838  ,, 


second  =  15.75 


^  r  nr 

4-1,  •  j  CO  15.75  T 

third  -  =  -  ■  =  1.095 

m  14 . 385 
po 


third  N 


mu 


1  +  1.095 
3.838 


=  0.5460 


From  figure  11  for  =  0.7 
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Third  =  0.8368 


APPENDIX  D 


COMPUTER  PROCEDURES  FOR  THE  GOVERNING  EQUATIONS 
WITH  SPHERICAL  GEOMETRY 


The  computer  procedures  described  for  spherical  geometry  in  this  appendix 
form  an  essential  part  of  the  so-called  "detailed"  analytical  model  but  can 
also  be  specialized  for  a  variety  of  less  detailed  models.  There  are  three 
computer  procedures . 


BASIC  COMPUTER  PROCEDURE 


The  first  computer  procedure  described  is  the  "basic”  procedure.  It^ 
cannot  be  regarded  as  a  design  since  it  contains  no  provision  for  automatic 
determination  of  crushable  casing  parameters  to  achieve  a  desired  acceleration 
and/or  a  desired  ratio  of  stroke  to  available  stroke.  Theprimary  purpose  of 
the  basic  procedure,  then,  is  to  check  the  adequacy  of  designs  which  have  been 
determined  by  other  means  (for  example,  designs  determined  by  the  simplified 
model  employing  figs.  5  through  12). 

With  symbols  defined  in  the  section  on  notation  or  in  parentheses,  the 
basic  procedure  contains  the  following  steps: 

1.  Input  and  print  the  following:  case  number,  WpQ,  Rp,  Uq,  £d>  ^mj 

^md»  Se’  ®L’  ®m  C^alue  of  g  for  equation  of  motion),  g^  (value  of  g  for 
acceleration  ratio)  ,  (value  of  e]c  for  equation  of  motion)  ,  e]ca  (value 

of  E]^  for  acceleration  ratio),  eks  (value  of  ek  fo^  stress  ratio),  Pckm^e 
(value  of  PckSe  equation  of  motion),  Pcka^e  (value  of  PckSe 

acceleration  ratio),  PcksSe  (value  of  PckSe  stress  ratio),  FpomC®) 

(value  of  Fpo(e)  for  equation  of  motion),  Fpoa(e)  (value  of  Fpo(e)  for 
acceleration  ratio) .  Also  input  and  print  trial  values  of  R  and  Oq. 

Finally,  input  and  print  whether  the  payload  is  considered  perfectly  bonded 

or  perfectly  unbonded  to  the  crushable  material,  whether  the  SEA  is  computed 
from  Oq  or  selected,  whether  the  material  is  considered  to  be  balsa- like  or 
honeycomb- like  (important  only  if  SEA  computed),  and  the  value  of  SEA  if 
selected.  Note  that  PckmSe’  PckaSe*  PcksSe  often  is  specified  as  PcmSe* 
to  be  calculated  in  step  2.  Note  also  that  Fpom(o)  PpoaC®)  specified 

as  a  constant,  usually  1.00,  or  as  the  integral  of  equation  (A45) . 

2.  Calculate  and  print  the  following  constants:  Ppge  (from  eq.  (A43) 
with  iHpo  =  WpQ/gg) ,  SEA  (from  eq.  (9)  if  calculated  for  balsa-like  or 

eq.  (10)  if  calculated  for  honeycomb -like) ,  PcmSe  (from  eq.  (8)  with 

factor  for  dimensions),  Wco  (from  Wco  =  (4/3)tt(r3  -  Rp3)pcmge),  +  Wpg, 

gm/i^-mdSe*  ga/’^mdSo  ^o^  (from  the  first  of  eqs .  (A33)),  %  (from  the  last  of 
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eqs.  CA33)  with  =  Wpo/g^  and  a  factor  for  dimensions),  and  Ppnge  (from 

eq.  (A35)  with  mpo  =  Wpo/g^) . 

3.  Generate  and  print  values  of  (from  eq.  (A36)  with  e]^xn  of 

step  1)  and  w^  (from  the  equation  of  motion,  eq.  (A37) ,  with  eq.  (A34)  and 

with  ejcm,  gm^  and  Pckm  of  step  1)  for  selected  values  of  z  (z  =  q/R)  . 
Terminate  the  integrations  at  the  lowest  positive  value  of  w^  when  z  is 
varying  by  increments  of  0.0001.  Print  this  value  of  z  separately  and  label 
it  without  penetration.  Calculate  L/R  according  to  equation  (5)  with 

^Pmax  ~  ^Pmax"^^^  print  it  separately,  and  label  it  L/R  without  penetration. 

4.  Generate  and  print  values  of  (from  eq.  (A36)  with  of 

step  1)  and  the  stress  ratio  (from  eq.  (A38)  with  eq.  (A34)  and  with 

e]^s  and  Pcks  of  step  1)  for  selected  values  of  z.  With  z  varying  by  incre¬ 
ments  of  0.0001,  terminate  the  calculations  at  z  =  1.00  or  at  the  lowest 
positive  value  of  1  -  (Opok/^o) >  whichever  occurs  at  a  lower  z  value.  Print 
the  termination  value  of  z  separately  and  label  it  zs.  Calculate 
Zpmax^^s^  plaint  it  separately,  and  label  it  Zp^^^/ Zg  without  penetration. 

5.  Generate  and  print  values  of  F^ia  (from  eq.  (A36)  for  eka  of 

step  1)  and  deceleration  ratio  np/n^d  (from  eq.  (A32)  with  eq.  (A34)  and  with 
^ka?  ga^  Pcka  from  step  1)  for  selected  values  of  z.  Terminate  the 
calculations  at  without  penetration. 

6.  If  the  payload  is  bonded  or  if  Zp^^^/z^  ^  1.00  without  penetration, 

terminate  the  program  at  step  5.  If  the  payload  is  unbonded  and  if 
Zp^ax^^s  ^  without  penetration,  proceed  into  the  penetration  phase.  For 

penetration,  calculate  the  constant  Kp  according  to  the  last  of  equations 
(A39)  and  calculate  four  initial  conditions  according  to  equations  (A40) . 

7.  Generate  and  print  values  of  y,  e,  dy/dx,  and  de/dx  for  selected 

values  of  x  (x  e  t/nnidge/%)  hy  integrating  the  two  simultaneous  ordinary 
differential  equations  (A41) ,  with  auxiliary  equations  (A42)  through  (A45) , 
with  initial  conditions  described  in  step  6,  with  Fciin(y)  determined  as  in 
step  3  when  (Rp/R)y  is  substituted  for  z,  with  Fponi(e)  being  the  integral 
of  equation  (A45)  or  a  constant  according  to  the  specification  in  step  1,  and 
with  ekm>  Pckm  being  the  values  of  step  1,  Terminate  the  integra¬ 

tion  when  dy/dx  is  zero  to  four  or  more  decimal  places  and  when  successive 
values  agree  to  four  or  more  signifcant  figures  for  the  worst  of  y,  e,  and 
de/dx.  ^The  corresponding  value  of  y  is  called  ymax^  and  the  corresponding 
values  of  e  and  de/dx  are  initial  conditions  for  the  next  phase  of  the 
problem. 

8.  With  the  initial  conditions  just  described,  generate  and  print  values 
of  e  and  de/dx  for  selected  values  of  x  by  integrating  the  first  of  equa¬ 
tions  (A41)  with  d^y/dx^  =  0.  Terminate  the  integration  when  de/dx  is 
zero  to  four  or  more  decimal  places  and  when  successive  values  of  e  agree 

to  four  or  more  significant  figures.  The  corresponding  value  of  e  is  called 
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emax-  Calculate  from  =  (Rp/R) (Xmax  +  ^kmemax) .  as  deducible  from 

equations  (A39)  with  Zp  =  qp/R.  Also  calculate  L/R  from  equation  (5),  and 
calculate  Print  zp^^^,  L/R,  and  zp^^^/zg  labeled  "with 

penetration, »» 

9.  Generate  and  print  values  of  the  deceleration  ratio  np/n^d 
selected  values  of  x  according  to  the  first  of  equations  (A41)  with  auxili¬ 
ary  equations  (A42) ,  (A43) ,  and  (A45) .  The  quantity  Fpoa(e)  is  the  integral 
of  equation  (A45)  or  a  constant  according  to  the  specification  in  step  1;  and 
£ka»  ga^  Pcka  values  of  step  1.  Step  9  terminates  the  basic 

procedure  when  penetration  is  present. 


SEARCH  FOR  OVERALL  RADIUS  R 


The  second  computer  procedure  described  is  a  design  procedure  in  which 
Oq  is  assumed  given  but  iterations  are  automatically  performed  to  determine 
R  for  a  desired  ratio  of  stroke  to  available  stroke.  The  desired  ratio  is 
unity  in  the  present  case,  that  is,  L/R  =  0;  but  the  available  stroke  leaves 
a  margin  of  safety  based  on  the  fictitious  compacting  strain  Cd* 


This  procedure  is  based  partially  on  a  modified  basic  procedure.  The 
modifications  include  the  specification  of  ^i^gg/n^j^d  Lhe  desired  value 

of  Up  ^/T^md^  determination  and  printout  of  (i*e.,  the  value  of 


at  n. 


‘P  "  "Pmax^  "“Pmax'“aes 

calculated  as  a  measure  of  the  acceleration  discrepancy 


and  B  (i.e.,  quantity  B  is 


The  next  step  is  to  iterate  the  modified  basic  procedure  just  described 
in  order  to  achieve  a  low  value  of  L/R.  As  a  start,  the  procedure  is  run  for 
two  values  of  R,  the  selection  being  based  on  the  design  charts,  figures  6 
through  9,  or  any  other  analytical  or  experimental  information  suggesting  low 
values  of  L/R.  The  results  are  labeled  (L/R) ^  for  Rj  and  (L/R) 2  for  R2. 

Then  the  iteration  is  based  on  successive  straight  lines  of  the  form: 

L/R  =  aR  +  b  (Dl) 

For  the  starting  values  (L/R) ^ ,  R^  and  (L/R)25  R2^  the  first  pair  of 
straight  lines  from  equation  (Dl)  is 


(L/R)i  =  a^^^Ri  + 

CL/R)2  =  a*^^^R2  +  b^^^ 

and  the  computer  determines  a^^^  and  b^^^  by  solving  the  two  simultaneous 
equations.  From  equation  (Dl)  for  L/R  =  0  (which  is  the  desired  value),  the 

computer  then  determines  R^^^  as 
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R 


CD  _ 


Then  is  run  through  the  modified  basic  procedure  to  determine  (L/R)  CD  . 

This  result  is  combined  with  the  two  starting  values,  and  the  two  having  the 
lowest  absolute  values  of  L/R  are  selected  for  a  new  pair  of  starting  values . 
The  process  is  repeated  until 


|l/R|  <  0.0005 

provided  that  each  new  value  of  1l/r|  is  lower  than  at  least  one  of  its  two 
starting  values . 

If  a  new  value  of  |l/r|  is  higher  than  either  of  its  starting  values, 
the  computing  machine  is  stopped,  and  new  starting  values  must  be  selected. 
This  did  not  happen  for  any  of  the  cases  reported  herein,  all  of  which  con¬ 
verged  rapidly  (requiring  an  average  of  0.17  min  of  execute  time  per  case 
without  penetration  and  4.55  min  for  the  detailed  model  with  penetration). 

In  contrast,  there  was  one  unrelated  case  that  did  not  converge.  For  this 
case,  however,  convergence  could  not  have  been  expected.  It  turned  out  that 
the  given  impact  velocity  was  too  high  for  a  feasible  energy  absorbing  design 
with  the  given  material,  and  hence  too  high  for  a  solution. 


SEARCH  FOR  OVERALL  RADIUS  R  AND  CRUSHING  STRESS  a 

o 


The  third  computer  procedure  used  herein  is  a  search  for  R  and  a  ,  which 
is  simply  an  extension  of  the  search  for  R  just  described  (except  that  pene¬ 
tration  is  not  included  since  eq.  (28)  determines  Og  for  penetration) .  In 
the  search  for  R  and  Og,  both  R  and  Og  are  varied  in  an  attempt  not  only  to 
make  L/R  =  0  but  also  to  make  6=0,  where  6=0  when  the  maximum  accelera¬ 
tion  load  factor  equals  the  desired  load  factor 

This  additional  requirement,  6=0,  makes  the  iterations  more  complicated 
than  before  although  the  modified  basic  procedure  to  be  iterated  is  the  same. 
This  time  the  starting  values  are  found  by  running  the  basic  procedure  for 
three  combinations  of  R  and  Og.  Again,  one  of  the  starting  values  can  be 
determined  by  the  design  charts  for  the  simplified  model,  figures  6  through  9, 
together  with  modifications  based  on  experience  for  other  models  and/or  other 
preliminary  information  (such  as  analyses  or  experiments  performed  on  similar 
configurations) .  Experience  indicates  that  the  other  two  starting  values 
should  be  small  deviations  from  the  first  in  which  the  higher  values  of  Og 
correspond  to  the  lower  values  of  R.  The  starting  values  are  labeled  (L/R)i 
and  6i  for  Rj  and  Ogi,  (L/R) 2  and  62  for  R2  and  002 >  and  (L/R) 3  and  63 
for  R3  and  ag3.  The  iterations  are  based  on  successive  pairs  of  planar 
surfaces  having  the  form 


L/R  =  aLaOo  +  appR  +  ap 


(D2) 


6  =  aggOg  +  ag^R  +  ag 


(D3) 
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For  the  three  sets  of  starting  values  just  labeled,  the  first  group  of 
three  planes  from  equation  (D2)  is 


D  1.  -  ^02  +  a  p  K2  +  a 


M  -  a^'^ 
R/3  '  La 


^La  ^03  +  Rs  +  a, 


and  the  computer  solves  the  three  simultaneous  equations  to  determine  a^^^ 

a*-^\  and  a^^^.  Similarly,  the  first  group  of  three  planes  from  equation  (D3) 
LR  ^  L 

31  =  aj^^aoi  +  a^^^^Ri  +  a^^^ 

32  =  +  ag^J^R2  +  a^'^ 

33  =  4o^Oo3  -  at^^Rs  +  a^'^ 

.  •  (1)  Cl)  j  (1) 

and  the  computer  determines  a^^  ,  a^j^  ,  and  a^  . 

Then  L/R  and  3  are  set  equal  to  the  desired  value  of  0  in  equations 
(D2)  and  (D3) ,  and  the  coefficients  just  determined  are  substituted  to  give 

^La  ""o  ^  \r  ^  \ 

^3a  ^3R  ^  3 

The  computer  then  solves  the  two  simultaneous  equations  and  determines 
(hopefully)  an  improved  pair  of  parameters  and  R^  ^ . 

The  parameters  and  R*-^^  are  introduced  into  the  modified  basic 

procedure  to  determine  the  corresponding  values  of  (L/R) and  3^^^.  This 
result  is  combined  with  the  three  sets  of  starting  values,  and  the  three  hav¬ 
ing  the  lowest  values  of  |l/r|  +  |3|  are  selected  for  new  starting  values. 
The  process  is  repeated  until 


|l/r|  <  0.0005 
|3|  <  0.0005 


with  the  requirement  that  each  new  value  of  |L/R|  + 
least  one  of  its  three  starting  values. 


is  lower  than  at 
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The  computing  machine  is  stopped,  analogously  to  the  R  search,  if  a 
new  value  of  |L/r|  +  |3|  is  higher  than  any  of  its  starting  values;  and  new 
starting  values  must  be  selected  (where  the  selection  can  often  be  facilitated 
by  a  plot  having  R  and  Oq  as  axes  with  values  of  |l/r|  +  |3|  indicated  by 
vectors  or  with  values  of  L/R  and  3  indicated  by  vectors) .  In  contrast  to 
the  R  search,  new  starting  values  had  to  be  selected  fairly  frequently  for 
the  search  for  R  and  Oq  (specifically,  for  roughly  one-fourth  of  the  cases), 
despite  the  fact  that  convergence  was  rapid  when  it  occurred  (0.7  min  of  exe¬ 
cute  time  per  case  for  the  detailed  model) . 


For  all  but  one  of  the  cases  requiring  a  new  set  of  starting  values, 
convergence  occurred  for  the  second  or  third  set.  The  exceptional  case  (not 
listed  in  table  1)  was  an  attempt  to  do  the  detailed  model  by  the  R  and  Oq 
search  so  as  to  match  an  np  of  2000  selected  for  the  simplified  model. 


For  the  landing  configuration  under  consideration 

different  values  of  Oq  indicated  that  an  n-n 

Pmax 


,  however,  R  searches  for 
as  high  as  2000  cannot  be 


attained  with  the  detailed  model.  Hence  the  iterative  method  would  have  been 


wrong  if  it  had  converged  to  an  of  2000  in  a  search  for  R  and  Oq. 
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APPENDIX  E 


EXACT  INTEGRATION  FOR  A  CLASS  OF  IMPACT  PROBLEMS  WITH 
VARIABLE  MASS  BUT  WITHOUT  PAYLOAD  PENETRATION 


The  class  of  impact  problems  to  be  considered  is  shown  in  figure  3(b). 

This  shows  a  general  vertically  symmetrical  landing  geometry  for  zero  shear 
resistance  and  uniform  compacting  strain  with  the  impact  to  be  described  before 
payload  penetration  occurs  (if  any) .  Equation  (A19)  defines  the  problem,  and 
all  assumptions  leading  to  that  equation  are  retained. 

One  additional  assumption  is  introduced,  namely,  a  uniform  and  isotropic 
crushing  stress  (o  =  o]^  =  constant)  .  Then  equation  (A19)  becomes 


i  d(U2)  =  g  dq 


>  ^  \o  -  fPck^c/"k^ 


(El) 


where  the  equation  has  been  multiplied  by  dq  and  where  Vc  has  been 
introduced  according  to  the  second  of  equations  (A14) . 

The  area  Acih  is  the  planar  area  of  crushable  material  that  is  flush 
with  the  landing  surface  in  figure  3(b),  and  Vq  is  the  volume  of  the  dotted 
region  in  figure  3(b).  Hence, 


^cih  dq  -  dV^  (E2) 

When  equation  (E2)  is  introduced  into  equation  (El),  each  term  in  the  latter 
equation  becomes  a  total  differential,  and  the  exact  integration  is 


1 

2 


U2 


gq 


^k^k 

^ck 


loge 


+  m 


CO 


(E3) 


where  C  is  a  constant  of  integration.  This  constant  can  be  evaluated  by 
noting  that  at  q  =  0,  U  =  Uq,  and  Vq  =  0 .  When  the  evaluated  C  is  intro¬ 
duced  in  equation  (E3) ,  the  result  is 


1. 

2 


-  U2)  = 


-gq  + 


^k^k 

Pck 


loge 


1 


1 


po 


m  ) 

CO 


where 

equations 


has  been  replaced  by  I  dh^ 

(A14)  .  •'o 


as  in  the  second  of 


(E4) 
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Equation  (E4)  gives  U  in  terms  of  q  and  can  be  simplified  if  the  gq 
term  is  neglected  (as  is  justified  for  most  impacts) .  Then  it  is  convenient 
to  reintroduce  and  let  it  be  the  only  variable  on  the  right-hand  side. 

Thus , 


1 

2 


(U„2  .  u2)  =  ^  log 


ck 


e  1  -  [PekV'kf”po  *  "co’l 


(E5) 


When  U  =  0,  Vo  has  reached  its  maximum  value,  .  Thus,  with  the  natural 

max 

log  converted  to  an  exponential : 


^c 

max 


♦  -cq) 

"ck 


(E6) 


It  is  implicit  in  the  derivation  of  equation  CA13)  that  Vq  +  Vci  = 
where  Vei  is  the  volume  of  m^  in  figure  3(b);  and  the  same  relationship 
applies  for  the  maximum  volumes.  Hence,  if  it  is  desired,  equation  (E6)  can 
be  written 


^  ^ cl 

^max  niax 


^o 


+  m 

CO 


1 


CE7) 


Equation  (E7)  agrees  exactly  with  equation  (1-4)  in  appendix  A  of 
reference  2  when  notational  differences  are  accounted  for.  The  development  in 
reference  2  includes  shear  resistance  implicitly  (by  showing  uplift  of  the 
compacted  region  in  a  figure)  but  counteracts  the  inclusion  of  shear  by  assum¬ 
ing  Aci  to  be  a  horizontal  planar  area  and  thereby  rendering  the  shear 
stresses  ineffective  in  energy  absorption.  Thus,  with  the  results  in  agree¬ 
ment  and  the  assumptions  reconciled,  the  development  of  equations  (E4)  through 
(E7)  becomes  a  partial  check  on  the  present  methods  as  well  as  a  set  of  exact 
results  for  a  special  case  of  variable  mass  in  the  absence  of  payload 
penetration. 
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TABLE  1.-  NUMERICAL  RESULTS  FOR  Uq  =  300  FT/ SEC,  =  0 , 7,  =  0 . 8 ,  AND  <  2000 
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(d)  Over  a  hypothetical  stroke  as  large  as  the  outer  radius,  there  is  no  force  developed  which  is  large  enough  to  cause  penetration. 


TABLE  1,-  NUMERICAL  RESULTS  FOR  Uq  =  300  FT/ SEC,  ed  =  0,7,  =  0.8,  AND  I  2000  -  Concluded 
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Stress 


Strain 


(a)  Stress-strain  curves. 


(b)  Crushing  test  with  plate. 


(c)  Crushing  test  with  plunger. 

Figure  1.-  Typical  stress-strain  curves  and  crushing  tests. 
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Payload 
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Figure  5.-  SEA  (for  =  0.8)  vs  Oq  for  selected  materials. 


70 


Figure  6.-  Stroke  design  chart  for  simplified  model 


Figure  8.-  Radius  ratio  design  chart  for  simplified  model  without  penetration 

when  Rn  is  given. 
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Figure  11.-  Mass  and  radius  ratio  design  chart  for  simplified  model  with  penetration  when  Rp  is  given 


Figure  12.-  Mass  design  chart  for  simplified  model  with  penetration  when  R- 
is  given;  also  simple  mass-volume  relation. 


Honeycomb- like  material 
Wpo  =  IOOIb  £^  =  0.8 

npmax  =  2000  -  0.7 

Uq  =  300  ft /sec 

Wco  SEA 


o  O  Simplified  model  1  ^  ^ 

□  A  Detailed  model  1  Penetration 


110.5  69.65  46.65  32.75 

Payload  packaging  density,  p.g. ,  Ib/ft^ 

f  P  c 

(a)  W(,Q  and  SEA 

Figure  13.-  Crushable  casing  properties  and  performance  as  functions  of  Rp 
and  ppge  for  honeycomb -like  material,  a  payload  weight  of  100  lb,  and  a 
payload  maximum  load  factor  of  2000. 
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Specific  energy  absorption  , 
SEA,  ftib/lb 


Overall  radius, 


Maximum  crushing  stress , 


value  of  payload  displacement,  qp 

-  >  - 

at  which  penetration  occurs  or  qs 

would  occur  if  unbonded 


See  part  (a)  for  design  conditions 
Pcm^e  ^Pmax 

o  bi  Simplified  model  with  penetration 

O  □  Simplified  model  1 

A  Detailed  model  ]  '^'^^out  penetration 
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Payload  packaging  density,  yO  g  ,  I 


‘IPmax/^s  PcmSe 


13  > 


Figure  13.-  Concluded. 


Balsa-like  material 
Wpo  =  450  lb 

SEA  =  24,000  ft  lb /lb 
Uo  =  300  ft /sec 


€n,  =  0.8 

6d  =0.7 

^  2000 


cTq  =  800  psi 

=  Ib/ft’ 

(Tg  =  1200  psi 

=5.76lb/ft^ 

rcm 


} 

} 


without  penetration 
with  penetration 


Undefined  transitions  for 


107.4  62.15  39.15  26.23 

Payload  packaging  density  ,  yOpgg,lb/ft^ 


(a)  Wco 

Figure  14.-  Crushable  casing  properties  and  performance  as  functions  of  Rp  and 
Ppge  for  balsa- like  material  and  a  payload  weight  of  450  lb. 
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See  part  (a)  for  design  conditions 
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(b)  R 

Figure  14.-  Continued 


Payload  maximum  load  factor, 


26.23 

g, ,  ib/f 


See  part  (a)  for  design  conditions 
_  Maximum  value  of  payload  displacement,  q 


Pmax 


Qo  at  which  penetration  occurs  or  would  occur  if  unbonded 


O  Simplified  model  1  yj^jf^out  penetration 
□  Detailed  model  J 

Simplified  model  1  penetration 
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Payload  packaging  density, 


Figure  14.-  Concluded. 
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